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I. INTRODUCTION 


The research grant NAG 1-500 entitled "Semiconductor Superlattice 
Photodetectors" was awarded to the University of Illinois at Urbana-Champaign 
by National Aeronautics and Space Administration-Langley Research Center 
on June 27, 1984. The grant was continued on July 1, 1985 and extended 
to December 31, 1986. Dr. Ivan Clark is the Technical Officer, and 
Mr. John F. Royall is the Grants Officer. 

This report is the final report. 

1. Period: 

June 1 , 1984 to December 31 , 1986 

2. Reporting Date: 

January 10, 1987. 

3. Technical Personnel: 

S. L. Chuang Assistant Professor of Electrical and Computer 
Engineering 

J. J. Coleman Professor of Electrical and Computer Engineering 

K. Hess Professor of Electrical and Computer Engineering and 

Research Professor of Coordinated Science Laboratory 
J. P. Leburton Assistant Professor of Electrical and Computer 
Engineering and Research Assistant Professor of 
Coordinated Science Laboratory 


Two research assistants 
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II. TECHNICAL PROGRESS 

Superlattice photodetectors have been investigated during the past 
two and one-half years. We have studied a few major physical processes 
in the quantum-well heterostructures related to the photon detection and 
electron conduction mechanisms, namely, the impact ionization of hot 
electrons, the tunneling-assisted process, the field effect on the wave 
functions and the energy levels of the electrons, and the optical absorption 
with and without the phonon assistance. The results are summarized 
below. The details are reported in the appendix where reprints and 
preprints of papers supported by this grant are included. 

1. Superlattice photomultiplier 

Tunneling-assisted impact ionization across the conduction-band-edge 
discontinuity of a quantum-well heterostructure is investigated and applied 
to a new superlattice structure (Appendix A). We consider multiquantum-well 
structures where the quantum-well regions are heavily doped and the 
undoped barrier regions are essentially insulating. Incident hot electrons 
due to the applied electric field perpendicular to the heterointerface 
interact with the two-dimensional electrons confined to the quantum wells 
through Coulomb force. The resultant electrons can either have enough 
energy to get out of the wells or tunnel through the triangular barriers. A 
new analytical approximation for the impact ionization rate is given which 
compares favorably with numerical results. The tunneling-assisted impact 
ionization rates and the ionization coefficients are calculated. It is 
shown that the tunneling effect reduces the ionization threshold and 
enhances the ionization rate significantly. Some experimental results of 
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this single-carrier type photomultiplier have been reported by Capasso's 
group at AT&T Bell Laboratories. 

We also carried out the variational calculations of subband eigenstates 
in an infinite quantum well with an applied electric field using Gram- 
Schmidt orthogonalized trial wave functions (Appendix B). The results 
agree very well with the exact numerical solutions even up to 1200 kV/cm. We 
also show that for increasing electric fields the energy of the ground 
state decreases, while that of higher subband states increases slightly 
up to 1000 kV/cm and then decreases for a well size of 100 A. 

We have also performed the exact numerical calculations for the energy 
level and the resonance width of quasibound states in a quantum well with 
an applied electric field (Quantum-well Stark Resonance) by solving the 
Schrodinger equation directly (Appendix C). This calculation gives both 
the resonance positions and widths for the complex eigenvalue E^, -if/2 of 
the system. Our theory also shows that the energy shifts of the ground 
states for the electrons and holes have the same behavior in high fields 
without any turnaround phenomenon, contrary to the results of Austin and 
Jaros . 

2 . Superlattice Photodetector Based on the Real Space Transfer Mechanism 

We have investigated the free carrier absorption process of a superlattice 
photodetector which makes use of the real space transfer mechanism. In 
particular, we have formulated and numerically computed the free carrier 
absorption coefficient for bulk GaAs in which we considered a second- 
order process involving both a photon and a phonon. Additionally, we 
carried out the computations for the free carrier absorption coefficient 
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in a quantum well for the above cited second-order process involving a 
photon and a phonon. 

We also carried out new calculations for the electric field dependence 
of the intersubband optical absorption within the conduction band of a 
quantum well (Appendix D). We show that for increasing electric field the 
absorption peak corresponding to the transition of states 1 -» 2 is shifted 
higher in energy and the peak amplitude is increased. These features are 
different from those of the exciton absorption. It is also found that 
the forbidden transition for states 1 + 3 when F - 0 is possible when F 
is nonzero. These results are significant for applications to infrared photo- 
detectors and infrared lasers making use of the intersubband transitions 


in the quantum wells. 
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Tunneling-assisted impact ionization across the conduction-band-edge discontinuity of 
quantum-well heterostractures is investigated and applied to a new superlattice structure. We 
consider multiquantum-well structures where the quantum-well regions are heavily doped and 
the undoped barrier regions are essentially insulating. Incident hot electrons due to the applied 
electric field perpendicular to the heterointerface interact with the two-dimensional electrons 
confined to the quantum wells through Coulomb force. The resultant electrons can either have 
enough energy to get out of the wells or tunnel through the triangular barriers. A new 
analytical approximation for the impact ionization rate is given which compares favorably with 
numerical results. The tunneling-assisted impact ionization rates and the ionization coefficients 
are calculated. It is shown that the tunneling effect reduces the ionization threshold and 
enhances the ionization rate significantly. 


I. INTRODUCTION 

Superlattice photodetectors have been of interest to 
many researchers because they involve new device physics 
and have the potential of minimizing the excess noise fac- 
tor. 1-4 The basic concept is based on the well-known fact that 
to minimize the excess noise due to the feedback process in 
an avalanche photodiode, one requires that the impact ioni- 
zation coefficients of the electrons and the holes differ great- 
ly, and the carrier type with the larger ionization coefficient 
initiates the impact ionization process. 5-7 Since the ioniza- 
tion coefficients, a a nd0, of the electrons and holes, respec- 
tively, can be very close, one tries to artificially enhance the 
a/0 ot 0 /a ratio using superlattice structures. 1-4 These in- 
clude a multiquantum-well structure, 1,2 a staircase photo- 
multiplier, 2,3 etc. One notes that most of these proposed de- 
vices deal with impact ionization across the energy band gap 
E t . We have studied recently the possibilities of a superlat- 
tice photomultiplier which makes use of the impact ioniza- 
tion across the conduction-band-edge discontinuity A £ c . 8 
Independent experimental work on this process has also 
been reported recently, 9 where both Al^In^As/ 
Gao 47 Ino .53 As and AlSb/GaSb superlattices are considered. 

The multiquantum-well device we consider is shown in 
Fig. 1, where the GaAs regions are doped with donors N D 
and the Al, Ga, _ x As regions are undoped. 8 Ohmic contacts 
are provided to the GaAs regions for replenishment of elec- 
trons in the wells once they are impact ionized by incident 
electrons which are heated up in the AiGaAs regions by the 
external electric field. A theoretical study of the impact ioni- 
zation in this structure has been done numerically in Ref. 8 
without considering the tunneli»$ effect. Since the impact 
ionization process is due to the eWciromconfined to the wells 
and virtually no holes are involved, tne electron ionization 
coefficient a will be much la^jzer than the hole ionization 
coefficient 0. In this paper, wrf also develop an anlytical ap- 
proximation for the impact ionization of electrons from the 


1 quantum wells and compare its results with those of our' 
previous numerical approach. We demonstrate very good 
agreement between the two methods. On this basis, tunnel- 1 
ing-assisted impact ionization across the band-edge discon- 
tinuity is calculated. Previous studies, 10,11 including the 
pioneering work of Keldysh and Kyuregyan, involve mainly 
tunneling-assisted impact ionization across the energy gap 
E g . We include Fowler-Norheim tunneling 12,13 by using a 
numerical approach involving Airy functions instead of us- 
ing the WKB method since the WKB method will not be 
appropriate when the energy of the electron is close to the 
barrier height (A£ e in this case). A comparison of the re- 
! suits using Airy functions and those for the WKB method is 
made. It is shown that tunneling reduces the threshold ener- 
j gy for impact ionization and enhances the ionization rate, 

* especially at the low-energy side. Numerical results for the 
electron ionization coefficient a are then presented with and 
I without the tunneling effect. The multiplication factor is also 
i shown versus the applied voltage. 



FIG. 1. Geometrical configuration of a multiquantum- well structure exhi- 
biting impact ionization across the conduction-band-edge discontinuity for 
! a solid-state photomultiplier. 
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In this section, we present (in Sec. II A) an analytical 
approximation for the impact ionization of electrons con- 
fined to quantum wells excluding tunneling effects, e com- 
pared the results with those obtained by using a direct nu- 
merical approach from Ref. 8. We then include the tunneling 
effect in the impact ionization process in Sec. II B. 

A. Analytical approximation for the average impact 
ionization rate 

The total transition rate per unit volume P tr between the 
incident hot electrons in the free states and the cold electrons 
confined to the quantum wells is given by 8 


has also been evaluated numerically in Ref. 8 where k 2 has 
both transverse component k 2 , and z-component 
±z(nir/L zn ). We know that k 2 , is limited by the doping 
concentration which determines the Fermi energy E F in the 
wells, 

j /tk*,)- 

I K% 1 = V(2 «•/*■)£, - (nir/L„ ? , (6b) 

I where n refers to the nth subband. The electron concentra- 
tion is 

n~N D = f f(E)g(E)dE = Y N„ , (7) 

Je qx 


1 ifO <k,,<Kir 

.0 otherwise 


'.-■pZZZZ 

V k ' k = hi hi 


X P^:/(k,)/(k 2 ) [ 1 -/(k; ) ] [ 1 -/( kj ) ] , (1) 


where 


2 IT, 


as -y-|<12|/f,|l , 2’)| 2 5(£ l +£,-£; —E'i). 


( 2 ) 


X exp 


\ 2 m*k a T,)i 


1 


(3) 


£*) 4 fie 2 V L J J ' 


S(k„k : 


.»-/ 

Jk u 4 - 


»«<* ih (2tr) 

S(E\ 4- E 2 — * E J — E f 2 ) 


[|k t — k;| 2 + 


k;-k, + kj 


( where 

i 

1 m* 


N. = 


The square of the matrix element 1 < 1 1'2'> ( 2 due to 

Coulomb interactions has been calculated in Ref. 8 where 
the initial state consists of electron 1 in free state above the 
well and electron 2 in bound state inside the well. The final 
state consists of both electron V and electron 2 ' in free states 
above the well. The average ionization rate ( 1/r) is given by 
Eq. (25) of Ref. 8: 


ME*) 

where n 0 is the concentration of the incident electrons, m * is 
the effective mass in the AlGaAs region, and 

r(E, 


(2 *) 2 L„ 

X /( k* ) [S( k,,k 2 : ) + S(k„k£; ) ] . (4) 

The summation is over all the occupied subbands in the wells 
and E 0z is the energy of the incident electron measured from 
the conduction-band edge of the barrier region. All param- 
eters are defined as in Ref. 8. An effective well width L zn is 
used to account for the finite barrier height. 8 The integral 
S(k,,k 2 ) due to the electron-electron interactions is defined 
as 

dk; 


(5) 

where the range of integration for k \ t is limited since both 
k \ t and k ^ must be >k ' h for electrons 1 ' and 2' to get out of 
wells if the tunneling effect is not included. The screening 
parameter q has been discussed in Ref. 8. The above integral 


ir&L. 


■& fl rin(l + e lEr ~ E °" )/k,T ) 


m' 


rrtfL, 


<E f -E Q h ) 


( 8 ) 


(for E f — E Qn >k B D is the electron concentration in the 
rtth subband and E Qn is the energy level of the /ith subband. 
Since ^k, is the momentum of the incident hot electrons 
which is usually much larger than the momentum of the 
electrons in the wells, and k, is essentially directed in the z 



r dkte . 

f it# 


\rl /: 0 J 

hr 1 

lm*k a T, J 




FIG. 2- A momentum space diagram illustrating the range of integrations 
for k; . Because of both energy and momentum conservation, the four wave 
vectors k t , k 2 , k; , and k 2 lie on a spherical surface in the momentum space 
assuming parabolic band structures. 
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Eq. (5) analytically noting that 
) 

— k,Mk 2 — k; ) 3 . ( 10 ) 

j and ignore the transverse components k J( in the integrand of and making a change of variable 
the integral S(k t ,k 2 ). Using the geometrical configuration s 

! shown in Fig. 2 in which k,, k 2 , k; , and k 2 lie on the same k = kj — k( , (11) 

j spherical surface 8 due to both energy and momentum con- j we obtain 


direction, k, — zk it we thus assume ' servation, we integrate 

5'(k t ,kj)»5(k |f lk Ji ), < 9 > &{E ' + E ' ~ E ' ~ E * 


S(zk u jk 11 ) 

V m* ^ (k ij ^21 — %k ) 

UirW/ [(fc„ — k lt )(k u -*,'„) + ? 2 ][(*i, -k u ) (k ; \,-k u ) + ? 2 ] 

= ' ( — \ jjji-Z**! if k • <k 

{tTffrq-JlK-k^ + q* ' h * 

0 otherwise 


_ it 


(12a) 

(12b) 

(12c) 


I In the first case, the condition k u <k ' h < (A:,, + k lc )/2 cor- 
responds to the integral over the shaded region in Fig. 2. If 
the threshold value k is <k 2l , the integral will be over the 
■ whole spherical surface 8 and results in the second Eq. ( 1 2b ) . 
j The threshold wave number k t ' h for electrons 1' and 2' is 
given by 

| • k; b =V2m?A T'/fi (13) 

which is larger than k u usually where the voltage drop 
| across the quantum-well region is assumed to be zero since it 
i is heavily doped and ohmic contacted. Here (or m*) is 
the average effective mass in the GaAs region with the non- 
! parabolic effect included. 8 Thus, the threshold energy of the 
! incident hot electrons for impact ionization £,, measured 
! from the conduction-band edge in the barrier region, is cal- 
i culated from 


j IONIZATION 
RATE {l/stc> 

)0 H ^ 



ENERGY Of INCIDENT ELECTRONS (tV) 


l = k\ x +*L -“^>2/c; h -k^ , (14) 

E, = {#/2m$2k ; b -k u ) 2 -&E c . (15) 

,! If one ignores k u = ir/L tX of the ground state, we have 

| £ ( =*3A£ C . Thus, the threshold energy £, is in general 

I =s 3 A£ c or less when k u is included. The integration over k 2f 

is easily done using the distribution function /(k 2 , ) in (6) 
and the fact that S( k , ,k 2 ) is independent of k 2 , . 

i i _ 7 wf£,y n" 1 

j t(E 0 ) . tfe 2 \ L ) 4irL„ 

; x K w 6 +, ( k *- , 61 - ,16) 

where E 0 = E ^ since the incident electron 1 is assumed to be 
traveling in the z direction only. The last integral over 
(or £ 0 ) in (3) to obtain the average ionization rate { 1/r) is 
calculated numerically by Simpson’s rule. We ignore k 2 , in \ 
the integrand of Eq. ( 5 ) which means that we assume that all 
electrons in the same subband in the quantum wells are ion- 
ized at the same rate independent of k 2 ,. Numerical results 
as shown in Fig. 3(a) indicate that this is an excellent ap- 
proximation which speeds up our computation by more than 
a factor of 100. The results using the direct numerical inte- 



FIG. 3. (a) The impact ionization rate ( \/r(E 0 ) J as a functiorrof the inci- 
dent electron energy E 0 measured from the conduction-band edge of the 
AlGaAs region using the analytical approximation (12) in Eq. (4). 
(circles — first subband, crosses second subband, dashed line = third 
subband, solid line ■= total. ) (b> l/r(£„) calculated usmg direct numerical 
integrations of multiple integrals as has been done in Fig. 4(b) of Ref. 8. 


grations 8 of multiple variables are shown in Fig. 3(b) for 
comparison. It is clear that our analytical approximations 
lead to very accurate results and eliminate the fluctuations 
due to numerical instabilities in the previous numerical inte- 
grations for all subbands as shown in Fig. 3(b), 


B. Tunneling-assisted impact ionization 

The impact ionization of electrons in the wells which 
subsequently tunnel through the triangular barrier is calcu- 
lated from 

f., -4:111 1 /fipWi )T(E;)A k.l/tk.) 

V ki kj 


x[i-/(k;)][i-/(k;)] f d7) : 

i 

where T(E \ ) and T(E' 2 ) are the transmission coefficients 
through a triangular barrier with a barrier height V 0 = A E c 
as shown in Fig. 4 at energies E \ and E 2 , respectively. The 
tunneling coefficient using the WKB method is given by 13 



where 

V Q =(K F /eF)(V 0 -E z ) f (19) ! 

K F = {2m'eF/#) m , (20) 1 

I 


! TRANSMISSION 
! COEFFICIENTS 



guiar barrier with a height LE e = 0.254 eV. The solid curve is using the 
numerical approach with Airy functions and their derivatives. The dashed 
curve is from the WKB method. 


where the effective mass difference between the well and the 
barrier regions has been ignored. If one includes the effective 
mass difference, the tunneling coefficient will be a function 
of both the longitudinal and the transverse energies, 13,14 E x 
and £ f . The effect is small, however, since the effective mass 
difference is small. A more accurate result for the tunneling 
coefficient without any asymptotic expansion can be ob- 
tained by use of Airy functions 13 : 


£(£,) = 


I 

*k z [Bi(77 0 ) - (Kf/k x )Aii'(7) 0 )] z + [Ai(77 0 ) + (K F /k z ] 2 ’ 


where 


I 

we have 


(21) 

! 

I 


k z = {2m*E z /fir) U2 (22) 

where tj 0 is defined in ( 1 9 ) , Ai ( rj 0 ) and Bi ( tj 0 ) are the Airy 
functions, and Ai' ( tj 0 ) and Bi' ( rj 0 ) are their derivatives with 
respect to the arguments. A numerical comparison of the 
two results using (18) and (21) is shown in Fig. 4 for a 
barrier height V 0 = A£ c « 0.254 eV. It is clear that the re- 
sults using the WKB method are good only if the condition 
|t 7 0 |>1 is valid. For energies E z near the barrier height, 
which are important for our calculation of the tunneling- 
assisted impact ionization, the WKB method is not accurate 
since the prefactor in (18) goes to zero, while the more pre- 
cise result shows a smooth transition. If one sets the prefac- 
tor equal to one as done by other researchers, 13,14 the dashed 
curve will shift to 1 at E = K 0 , and will deviate more from the 
solid curve, which is obtained from the exact numerical ap- 
proach, than the dashed line. Thus, we conclude that (21 ) 
should be used to incorporate the effect of the tunneling- 
assisted impact ionization. 

Since the tunneling coefficient is a function of E z only, 


S(z*„,k 2 ) =tS(zk lgt zk u ) 

m* p' 
Kv 2 # Jk „ 


dk\ z 


X 


nE\ z )T{E^) 


(23) 


[(*!, -*!«)(*!, -K ) +<r] 2 

where T{E\ Z ) and r(£^ ) are obtained from (21 ) and the 
range of integration is from k ^ to k lzf i.e., over the whole 
spherical surface in Fig. 2, instead of from k[ h to 
k\ z + k^ — k ; h since the tunneling effect is taken into ac- 
count in T(E\ Z ) and 7(£^). The previous calculations 
without tunneling correspond to 


£(£;,) 


1 if£;,>A£ c 

0 if E \ 2 < A£ c 9 


(24) 


and a similar expression holds for T{E , Zi ). Thus, the results 
of impact ionization without tunneling will have threshold 
energies decided by ( 15), which are higher than that of the 
tunneling-assisted impact ionization since the transmission 



coefficient differs significantly from zero even if E \ z and 

E^<&E C . 


III. NUMERICAL RESULTS 

Various effects of multiple subbands, the quantum-well 
width L x , the band-edge discontinuity LE e1 and the doping 
concentration have been discussed in Ref. 8. We have used 
the 60-40% rule for the conduction- and valance-band-edge 
discontinuities. 15 The validity of our new analytical approxi- 
mation in Eq. ( 12) has been presented in Fig. 3(a) and com- 
pared with the direct numerical integrations shown in Fig. 
3(b) using Eq. (4). The results of the analytical formula in 
Fig. 3(a) do not show the fluctuations of the multiple nu- 
merical integrations in Fig. 3(b). 

In Fig. 5, we show the ionization rate as a function of the 
incident electron energy E 0 measured from the conduction- 
band edge of the barrier region at the heterointerface with 
and without tunneling. The doping concentration N D is 
1 X 10 17 cm“ 3 and only one subband is filled in this case. It 
can be seen that the threshold energy for impact ionization is 
reduced significantly when tunneling is included. The ioni- 
zation rate is also enhanced especially at the low-energy end. 
When the energy of the incident electron is very high, the 
tunnelings of the final electrons V and 2' are not important 
which is expected, since they all have enough energy to get 
i out of wells without tunneling assistance. In Fig. 6, we calcu- 
i late the ionization coefficient a for a doping concentration 
■V D = 7 x 10 18 cm' 3 and assume the saturation velocity of 
the electrons v s = 1 x 10 7 cm/s. The ionization coefficient is 
obtained using a = ( \/r)/v s . The electron temperature T, 
versus the electric field F is obtained from the average energy 
: by a Monte Carlo simulation and is roughly, e.g., T, = 4630 
j K at F = 100 kV/cm; T ( = 3860 K at F= 50 kV/cm; and 
J T t = 2300 K at F= 10 kV/cm. The multiplication factor 
1 M n for electrons as a function of the applied voltage V for a 
: multiple-quantum- well structure of 50 periods is shown in 
Fig. 7. We assume a well size L z — 200 A, and a period 


IONIZATION 
COEFFICIENT (I/cm) 


of^ al page is 
F *** QVAL,n 



0 


X 



t 


X 


Tunneling Assisted 
Impact Ionization 


Without Tunneling 


5 



1.0 2.0 3 0 4 0 

I/ELECTRIC FIELD U l0' 5 cm/V) 


FIG. 6. I ne electron ionization coefficients a with and without tunneling 
) assistance vs the reciprocal of the electric held are shown. 


j L = 600 A. At this high-doping concentration, there are 
I four occupied subbands. The average effective mass mj 
! = 0.09 1 m 0 for the GaAs region is used to take care of the 

I effect of the nonparabolicity. 8 We have/n * = 0.083 m 0 in the 
I AljGaj^As region assuming the Al mole fraction 
x = 0.25. The conduction-band-edge discontinuity is A£ c 
= 0.187 eV. 

All the calculations are done assuming T = 77 K since 
the confinement of the carriers by the quantum wells will be 
j -better at low temperature. This dynamic storage of carriers 9 
| in the quantum wells is very crucial in the design of this 
! device. According to Eqs. (4) and ( 12), the impact ioniza- 
\ tion can also be written as 

1 _ y N* 4ire 4 (L zn V 

r(£ 0 ) ~ $ 2 ^ V L J 

I x K k '- i 6 +s (‘"- 1 i:)]’ <25) 
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FIG. 5. The impact ionization rates [ l/r(£ 0 ) ] as functions of the incident 
electron energy with tunneling-assisted process (circles) and without tun- 
neling (crosses) are compared. 
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FIG. 7. The multiplication factor for electrons vs the applied voltage in 3 
multiquantum-well structure with 50 periods is illustrated above. 


i.e., the impact ionization rate due to electrons in each sub- 
band of the quantum well is proportional to the electron 
concentration in that subband. This is true provided that 
S(k,,k 2 ) ^5(k,4A'2j )» » the transverse component k 2 , is 
negligible in the 5 integrals. To increase the ionization rate, 
one either increases AlE c or increases L z to have more elec- 
trons N d stored in the wells. However, larger A£ c will in- 
crease the threshold energy, and larger L z will decrease the 
energy of the confined electrons. Both tend to decrease the 
impact ionization rate. Thus, one needs to choose optimum 
parameters for A E ct L xt and doping concentration N D for 
best performance of the device. The multiplication factor for 
electrons is calculated from 

M„ = (l+a£)*, (26) 

where the 1 accounts for the primary electrons, aL accounts 
for the secondary electrons and N is the number of periods. If 
the total length of the avalanche region is W , we have 
W = NL y and in the limit A T -> oo 

M n = lim [ \+a{\V/N)) s = e aW , (27) 

which is the same as that of an ideal avalanche photodiode 16 
with an avalanche region W. However, since L is finite, we 
use ( 26 ) instead of ( 27 ) . The multiplication factor increases 
steadily or almost exponentially for this superlattice photo- 
multiplier instead of increasing sharply at a certain reverse- 
biased voltage which would be the case for a two-carrier- 
type avalanche photodiode. 16 

IV. CONCLUSIONS 

Tunneling-assisted impact ionization across the con- 
duction-band-edge discontinuity has been investigated theo- 
retically. An excellent analytical approximation has been 


found for electron -electron interactions between the inci- 
dent free electrons and the bound two-dimensional electrons 
in the quantum wells. It is shown that the tunneling-assisted 
process lowers the threshold energy and enhances the ioniza- 
tion rate. The tunneling may increase the ionization coeffi- 
cient by a factor of 2 or 3 depending on the electric field 
strength and the size of the quantum wells. 
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We present variational calculations of subband eigenstates in an infinite quantum well with an 
applied electric field using Gram-Schmidt orthogonalized trial wave functions. The results 
agree very well with the exact numerical solutions even up to 1200 kV/cm. We also show that 
for increasing electric fields the energy of the ground state decreases, while that of higher 
subband states increases slightly up to 1000 kV/cm and then decreases for a well size of 100 A. 


Electronic as well as optical properties of quantum wells 
subject to external electric fields have received much atten- 
tion. 1-9 These areas are interesting both from a fundamental 
and a practical point of view. The possible device applica- 
tions include the use of quantum confined Stark effect 2 in 
optical modulators* and optical switching devices. 9 As yet, 
most of the theoretical works has been confined to the calcu- 
lations of the ground state. Very recently, Matsuura and Ka- 
mizato" reported an exact numerical calculation of subbands 
in an infinite well and concluded that the higher subbands 
behave very differently from that for the ground state when 
the electric field strength is increased. Exact solutions em- 
ploying two independent Airy functions 2,5 are sometimes 
too complicated to use in real problems. The variational ap- 
proach has the advantage of providing analytical expres- 
sions for the eigenstate energies and the wave functions, and 
numerical results with reasonable accuracy can be obtained. 
Analytic form of the trial wave function for the ground states 
has been known/ 2 but no useful forms of the trial solutions 
for the higher subbands which yield accurate results com- 
pared with the exact ones are given yet. Recently, a trial 
wave function of the form /„(z)exp( — fi c 2 ) for the nth 
subband has been suggested/ where f „ (z) is the zero-field 
nth quantum well bound state wave function and 0 e is a 
variational parameter. However, it is pointed out that the 
solutions may yield very different results from the exact ones 
(for n > 1 ) because these trial functions are not orthogonal 
to each other. 7 That observation is also confirmed in this 
paper and the numerical results are illustrated. Thus, it is 
important to find an orthogonalized set of trial wave func- 
tions for the variational calculations of subband energies and 
wave functions. 

We report in this letter variational calculations on sub- 
band states in an external uniform electric field based on the 
infinite well model. We find analytic forms of orthogona- 
lized trial wave functions by the Green-Schmidt orthogona- 
lization procedure. Our calculations agree very well with 
exact numerical results up to 1200 kV/cm with an error less 
than 9% for L = 100 A. For comparison, we also show the 
numerical results for which the trial wave functions are not 
orthogonalized. 

It is well know n that the variational method can also be 


used to obtain one of the higher energy levels if the function 
is orthogonal to the eigenfunctions of all the lower states. 10 
The most common method of obtaining an orthogonal set of 
functions is the Gram-Schmidt orthogonalization proce- 
dure, 11 which is the construction of an orthonormal set 
{^,,^ 2 ,...} from a finite or an infinite independent set 
which is not necessarily orthonormal. 

Let us consider an electron with charge — |e| and effec- 
tive mass m*, in an infinite quantum well width L in the 
presence of a constant electric field F along the positive di- 
rection of the well z. We choose the origin to be at the center 
of the well. The Hamiltonian of the system in the effective 
mass approximation is given by 1,2 

H = H 0 +\e\Fz t (1) 

where H 0 is the unperturbed Hamiltonian whose eigenvalues 
are given by 

£<°> = (*V/2m*£ 2 )n 2 n = 1,2,... (2) 


For our specific problem, we chose the nth vector u n to be 
nir i 


M„(z) = sin 


X exp 


rK)l 


i ■ L 
l z l <y > 


(3) 


which is not an orthogonal set, where is the nth variation- 
al parameter. One can easily see that for n = 1 , u , is the trial 
solution introduced by Bastard et aL ] We also define the 
inner product between two functions / and g , if\g), by 

</|g>=f f*(z)g(z)dz, (4) 

J - L/2 

where the superscript • denotes a complex conjugate. 

The procedure we use is as follows: 

n = 1, (i) Let tA, = w, . (5a) 


(ii) Minimize £,(/?,) = <tf ,| H ,(5b) 

and find /?,. 

(iii ) <A, = !/•,/< t/' f | i/'i) 1/2 . (5c) 

n — 2, (i) Let tlt 2 = u z — {u 2 \d>\)d>\ . (6a) 


(ii) Minimize £ : (/? : ) = (t/'/// ,(6b) 
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and find /? 2 . 

(»i) = \ M<)MW /2 - 

t - I 

n = l, (i) Let ip, = u, - 


(ii) Minimize £,(£,) = <#,|//|* f )/<*,|*,) , (7b) 
(6c) and find)?/. 


(7a) 

I 


(iii) <p, = 


1/2 


(7c) 


This procedure requires mainly three inner products which are in analytical forms: 

, , v L, 0 . 0 J\-(-\)- m txp[-tf n +0 m )] 1 -( -!)" — exp[-(/?„-f/? m) n 

{u - \u m ) - T (ft +e m ^ ( * + „) v + <a + )= ; ’ 

±«r([« - »• ^->1 


(it - m)V + (£. + £„,) ! 


[,«*«>(»= 6 -) 1 !-(-!)— «p[- (ft +/>.>] ) 

-^rj-t»e. (» + i")j ( „ +m) v + <*+«.)■ j- 


and 


Wm> 2 |e|/X2 (« + m)V + ()? ) ,+iS m ) J 


[( 


1 + 


Pn + Pm 


2 (fi„ +P m Y 


7) 


(n + m) 2 ir + ( P„ +P„) 

X(1 - ( - 1)" + m exp[ — (P„ + P m )]} + (P„ +P m )( - l)" + m exp[ - (0„ +0 m )]] 

P„+P„ 2(P n 4- P m ) 2 


- — |e|FL 2 r-r-! 

2 (n — m) 2 ir + (P„+P m ) 


4( ,+ 


7) 


2 (n — m) : 7r + (/?„ +/?, 

x{l — ( — 1 )" “ m exp[ — (P„ + 0 m )]} + (0„ +0 m )(-I)""exp[ - Off. +&,)]]. 


( 8 ) 


(9) 


( 10 ) 


The minimization of an analytical function and finding 
/? m,n can be directly done in the computer by calling some 
subroutines, e.g., using International Mathematical Statisti- 
cal Libraries (IMSL) subroutines. 

For the expectation value of the ground-state energy 
£,(/?,), we have 

E\(P \ ) = |<£,) 


-MO + t 9 


__i.coth/?^j 


+ y {w,+j?ir 

where the parameter F is defined by 
£ = \e\FL/E\°\ 

which turns out to be the normalized electric field. 7 
In the low field limit, F< 1, 




and 


(ID 


( 12 ) 


(13) 


a ,£,=£, (/? r)-£T= ~ 7(7 ~ A)' ? {r ~ 4 


In the high field limit, /> 1 
0 7" = (JttF)'' 5 


(14) 

(15) 


and 


Equations (11 )— ( 16) have already been obtained by 
Bastard et al . 1 except they missed the factor of 1 /2 in front of 
0 in Eq. ( 7 ) and the expressions for /? ™ n of Ref. 1 in the low 
and high field limits are not correct. This does not affect their 
final results for which are correct. Calculations of high- 
er energy levels are straightforward following the proce- 
dures in Eqs. ( 5 ) - ( 7 ) . The results of the normalized energy 
E n = E n /E j 0) for the first three states are plotted versus the 
normalized electric field F in Fig. 1. The plot in terms of j 
these two normalized parameters is universal which can be j 
seen from the original Schrodinger wave equation. 7 It is ' 
readily seen that the shift of the subband energy due to the ! 
electric field is different between the n = 1 state and the j 
higher energy states. For increasing electric fields, the 
ground state shows a large negative shift, while higher states 
have small positive shifts for fields up to £=20 and then 
negative shifts. We have also plotted the subband energies 
obtained from the exact numerical solution. We find that | 
both methods of calculation give very close results, even up 
to F = 25. The parameters W'e use are m* — 0.065m f , and 
L — 100 A where m 0 denotes the free-electron mass. F = 10 
corresponds to 578.5 kV/cm for this L. The results of calcu- 
lations employing unorthogonalized trial functions^ defined 
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?--HFL/Ei 01 {NORMALIZED FIELD) 


0 5785 MV.O 1735.5 

FUV/cm) FIELD (for L*KX)A) 

FIG. 1. Normalized subband energies (£, = £„ /£j 0) ) are plotted vs the 
| normalized electric field ( F — \z\FL /£ S 0) ) . ( 1 ) The exact solution — ( sol- 

id lines), (2) the variational method with Gram-Schmidt orthogonaliza- 

tion procedure (dashed lines), (3) the variational method (see Ref. 5) 

| without orthogonalization (dotted lines). 

in Eq. (3) arc plotted in Fig. 1 (dotted lines). One readily 
j notices that these solutions begin deviating significantly 
I from the exact ones at aroundF = 1.5 (F= 86.8 kV/cm for 
L— 100 A). In Fig. 2 we plot the square of the wave func- 
| tions |0| 2 for the first three states when the electric field is 
F — 20. 1 0| 2 is normalized to L. We see in Fig. 2 that the 
ground-state wave function is shifted to z < 0 region signifi- 
cantly, While higher subband wave functions still have a 
large amplitude in the z > 0 region even at F = 20. 

In conclusion, we have derived orthogonalized trial 
i wave functions which yield very close results to those of the 
exact numerical solutions. It is shown that for decreasing 
electric fields the energy of the lowest subband decreases, 
while that of the higher subbands increases slightly up to 
1 20 (F> 1000 kV/cm for I = 100 A). 


ii*2 



i/L 

FIG. 2. Magnitudes of the wave functions, L (z)| 2 , are plotted vs the 
normalized distance z/L for the first three suhbands. The dotted lines are 
for the zero electric field. The solid lines are the exact solutions and the 
dashed lines are those for the variational method with the Gram-Schmidt 
orthogonalization procedure for a quantum well with an applied electric 
field (£« 20). 
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Exact calculations of quasibound states of an isolated quantum well with uniform electric field: 

Quantum-well Stark resonance 
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We present universal plots from exact numerical calculations for the energy level and the reso- 
nance width of quasibound states in a quantum well with an applied electric field (quantum-well 
Stark resonance) by solving the Sch rod infer equation directly. This calculation gives both the 
resonance positions and widths for the complex eigenvalue £«— /I72 of the system. Our theory 
also shows that the energy shifts of the ground states for the electrons and boles have the same 
behaviors in high fields without any turnaround phenomenon, contrary to the results of Austin and 
Jarcs. 


Electronic and optical properties of quantum wells with 
applied external electric Helds are of increasing interest 
Studies of these areas are important both from a funda- 
mental and a practical point of view. Optical modulators 1 
and optical switching devices 2 based on the quantum con- 
fined Stark effect have been suggested. Possible device ap- 
plications of the field-induced tunneling in quantum-well 
and quantum-barrier heterostructures include high-speed 
resonant tunneling devices.** 4 
More recent theoretical studies 7 ** of the effects of exter- 
nal electric fields on the quantum-well systems have 
predicted both the field-induced level shifts and the field 
dependence of the carrier lifetime. In this paper, we report 
exact numerical calculations on quasibound states of a 
quantum well in an external electric field (quantum-well 
Stark resonance) by solving the Schrddinger equation for 
Stark resonance directly. It is found that the previous re- 
sults based on phase-shift analysis 7 ** and the stabilization 
method* agree very well with our results over a wide range 


Potentioi Ciw*ey 



FIG. 1. Potential-energy profile Viz) for i single quantum 
well with depth Vq and width L subject to an external electric 
field F. 

a 


of the electric field. At an extremely high electric field, 
there is no turnaround behavior in the energy shift for both 
the electrons and the holes, contrary to the results in Ref. 
7, where no explanation can be provided for that 
phenomenon. We believe that our direct numerical ap- 
proach is very reliable even at a very high electric field, 
while the results using the phase-shift analysis may have 
drawbacks in the high-field limit Our approach has an 
advantage over the previous results 7 ** in that both the 
Stark resonance position (quasibound-state level) and the 
width can be obtained from the single complex energy 
eigenvalue of the quantum-well Stark resonance problem. 
The disadvantage is that numerical subroutines of the Airy 
functions with complex arguments are required. 

Consider an electron with charge - | * | and effective 
mass m*. in a finite quantum well of width L and depth Po 
in the presence of a constant electric field F along the posi- 
tive direction of the well z (Fig. 1). We choose the origin 
to be at the center of the well The Schrodinger equation 
of the system in the effective-mass approximation is given 
by 7 ** 

—£-&*(» >+M FrwfzWtfr). M *1/2 . 

2m dz 2 

a* d* * 

■ •- ) + (Ko+ 1*1 Fz )*<z ) -£y(z ), 

2 m* dz * 

\z\>L/2 . 

Since the potential energy term in Eq. (1) tends to — «• as 
z goes to — «•, the system does not, strictly speaking, have 
true bound states. 7,1 * In other words, the particle initially 
confined in a well can always lower its potential energy by 
tunneling out of the well when the field is not zero. It may 
happen, however, that the tunneling probability is very 
small In such a case, we can regard the system as having 
quasibound states, in which the particles move “inside the 
well” for a considerable period of time and leave through 
tunneling only when a fairly long time interval r has 
elapsed. In discussing the quasibound states, we may use 
the following formal method. Instead of considering the 
solutions of the Schrodinger equation with a boundary 
condition requiring the finiteness of the wave function at 
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infinity, we shall look for solutions which represent outgo* 
ing waves at infinity; 10 this implies that the particle finally 
leaves the well by tunneling. Since such a boundary condi- 
tion is complex, we cannot assert that the eigenvalues (en- 
ergy) must be real By solving the Schrodinger equation, 
we obtain a set of complex eigenvalues, which we write in 
the form 

£-£o~iT/2, (2) 

where T is found to be positive. £o and I* correspond to 
the quasibound-state energy level and the resonance width, 
respectively. The tunneling probability per unit time is de- 
fined by 

«-r Zh . (3) 

The solutions to Eq. (1) with the outgoing-wave condition 
are linear combinations of two independent Airy func- 


11 





«ilBi(ij 2 )+i Ai(i72)l, z < —L/2, 


yr(z) “■ 

aoAi(i 7 i)+b 0 Bi(i 7 i), |z|iSL/2, 
azAKife), z > -L/2 , 

(4) 

7i““ 

2m* 

1/3 

(£— |*|£z), 

(5a) 

.UhF)* 


2m* 

1/3 

(£—Vo~ |e |£z) . 

(5b) 

( ehF ) 2 


The wave function for z < —L/2 represents an electron 
traveling to z “ — «• after tunneling. The complex energy 
£ can be found by solving the secular equation obtained by 
matching the value of y and its first derivative at the 
points, z “ ±1/2. The resulting determ inantal equation 
is 


( 6 ) 


det 


Ai(jii + ) Bi(in + ) — Ai (q^) 
Ai'(»j| + ) Bi'(iji + ) — Ai'(ijj') 
Ai(rji — ) BiGjf) 0 

Ai'(qf) Bi'(jjf ) 0 


0 

0 

— lBi(i»D+/ Ai(ijj* )] 

- lBi'(ifc“ ) +/ Ai'(ijj” )] 


-0 


where iff* and rj? are the values of r/i and 172 evaluated at 
z “L/2 and —L/2, respectively. If we introduce a new 
parameter £ <0> defined by 

\2 




(7) 


(which happens to be the ground-state energy of an infi- 
nite quantum well with width L ), and define the normal- 
ized energy £“£/£ <0> , the normalized electric Held 
£“ |e |£L/£ (0> , and the normalized well depth K 0 
-K„ /£<«, we may express 7 * and 72 * by these three 
normalized quantities: £, F, and V<y 


m 


■ - Of* 


F 1 


1/3 

(ET+F) 


x 

F l 


1/3 

- (£ -K 0 T +£>- - 


(8a) 

(8b> 


TABLE I. Comparison of the numerical results for Et—iV/2 
using the exact numerical method of this paper, the phase-shift 
analysis (Refs. 7 and 9), and the stabilization method (Ref. S). 


F(kV/cm) 


This paper 
(eV) 

Phase-shift 
analysis (eV) 

Stabilization 
method (eV) 

75 

£e 

0.025167 

0.025 167 

0.025 167 


r 

1.86x10"* 

1.9x10"* 

8.6x10"* 

100 

£0 

0.0242107 

0.0242105 

0.0242106 


r 

3.64x10"* 

3.6x10"* 

4.1x10"* 

ISO 

£0 

0.0213716 

0.0213816 

0.021 170 


r 

6.41x10"* 

6.4x10"* 

6.5x10"* 


I 3 

This means that the solution of £ from Eq. (6) is universal 
and can be used for both electrons and holes with the re- 
placement of the parameter £ <0> with their corresponding 
effective masses. 1 ' (Here the effective masses inside and 
outside the well are assumed to be equal.) The normalized 
energy £ can be expressed in terms of only two normalized 
parameters, Vq and F. Thus it is clear that both electrons 
and holes should have the same behaviors in their energy 
shift and the resonance width. To obtain the results of £ 0 



FIG. 2. Comparison of the ground-state energy of the varia- 
tional calculation (Refs. 13 and 15) for infinite-well with ap- 
propriate effective-well width (dashed line) and the real part of 
the energy eigenvalue £a from exact calculation (solid line) of 
this paper. 
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FIG. 3. The real part of the normalized energy £i«£t/£ #l 
for various normalized well depth* is plotted vs the 

normalized electric Held / — |# | FL/E m 


and T for holes, one need only multiply £ by £o using the 
effective mass of the hole. We have solved Eq. (6) numer- 
ically to the desired accuracy using the series and asymp- 
totic expansions of the Airy functions with complex argu- 
ments. 11 To check the validity of our approach, we com- 
pared our results with those of the previous methods 7,1 in 
Table I. The values of Vo, L , and m* for the heavy holes 
used in the calculations are, respectively, 

Ko-100 meV, L -37 A, m* -0.45m 0 , (9) 


where mo is the free -electron mass. It is readily seen that 
our results agree very well with those of the phase-shift 
analysis 7 '* and the stabilization method.* 

In Fig. 2, the real part of the energy £o (resonance 
position- solid line) for the ground-state energy with the 
values of Vo, L, and m* for electrons given by Ko— 340 
meV, L — 100 A, and m* — 0.0665mo is compared with the 
results of infinite-well variational calculations 11 * 15 (dashed 
lines), where we have used an effective well width 
£«rr— 126.5 A, chosen to give the same £o at zero field for 
the variational calculations. It can be easily noticed that 
both calculations gave very similar results even up to 
2x10* V/cm. However, the variational calculation for the 
infinite-well model cannot give the resonance width since 
no tunneling exists for the infinite welL The results of the 
normalized resonance energy £o—£</£ <0) for various V 0 
are plotted versus F in Fig. 3. In contrast to the previous 
results 7 which are still controversial, the resonance position 
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FIG. 4. The normalized resonance width f — r/£ w> for vari- 
ous » plotted *« normalized eleeric field F — ! t j FL!E m 


is found to be in the well even at very high field. The 
behaviors of the resonant position are the same for both 
electrons and holes with proper £ (0> used together with 
Fig. 3 as discussed before. Thus the turnaround behavior 
for the holes and electrons in the energy shift shown in 
Ref. 7 is probably a drawback of that method itself. Using 
the same numerical values for holes as those in Ref. 7. 
£-30 A, Fo“ 70 meV, m* — 0.45mo, we obtain £ (0 ' 
—92.26 meV, Ko— 0.76. We do not have any turnaround 
behavior even up to F — 10, or the electric field F —3075 
kV/cm, which covers a much wider range of electric field 
than that of Ref. 7. In Fig. 4, we plot the normalized reso- 
nance width f — T/£ (0> for various Vo vs F. Since the life- 
time r is defined by r— h/T, the results plotted in Fig. 4 
predict a rapid decrease of the carrier lifetime with in- 
creasing applied field by field enhanced tunneling. 

In conclusion, we have solved the Schrodinger equation 
for a quantum well with uniform electric field directly. 
Complex eigenvalues for quantum-well Stark resonance 
are obtained. Our approach has an advantage over previ- 
ous analyses 7 '* in that both the resonance position and 
width can be obtained from a single complex eigenvalue of 
the problem. . - 

This work was partially supported by the Air Force, 
Contract No. F33613-84-K-1S57 and the NASA Grant 
No. NAG 1-500. One of the authors (D.A.) would like to 
acknowledge the support of GTE. 
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Abstract 

We present new results for the electric field dependence of the intra- 
subband optical absorption within the conduction band of a quantum well. 

We show that for increasing electric field the absorption peak correspond- 
ing to the transition of states 1 2 is shifted higher in energy and the 

peak amplitude is increased. These features are different from those of 
the exiton absorption. It Is also found that the forbidden transition for 
states 1+3 when F * 0 is possible when F is nonzero. 


2 


Quantum confinement of carriers in a semiconductor quantum well leads 
to the formation of discrete energy levels and the drastic change of opti- 
cal absorption spectra* * The interband absorptions near the band gap have 

been extensively studied, and it has been shown that their absorption and 

2 3 

luminescence spectra are dominated by excitonic effects. 9 Some more 

recent studies have concentrated on the electric field dependence of energy 
4-7 

levels and band edge optical absorption including the exciton 
8—10 

effect. Very recently, the experimental studies of the intrasubband 

absorption within the conduction band of a GaAs quantum well without the 

applied electric field have been reported**. A very large dipole strength 

and a narrow band width were observed. In this paper, we present the 

theoretical calculations for the electric field dependence of the optical 

absorption between the discrete subbands within the conduction band of a 

quantum well based on the infinite potential barrier model. One of the 

reasons for increased interest in this area is the possibility of practical 

12 

device application. For example, in 1970, Kazarinov and Suris proposed a 
new type of infrared laser amplifier using the intrasubband transition and 
resonant tunneling. Far-infrared photodetector with high wavelength selec- 
tivity based on the intrasubband absorption and the sequential resonant 

tunneling has also been suggested. *^ 

The Hamiltonian of the system (for a single quantum well) subject to a 

uniform electric field perpendicular to the quantum well (the z-direction) 

in the presence of an optical radiation (Fig. 1) is written as 


H 


H + H" 
o op 


9 


( 1 ) 


3 


where H q is the unperturbed Hamiltonian for an electron in the quantum well 

in the presence of perpendicular electric field, and the interaction 

>* 14 

Hamiltonian H is given by 
op 


H 


op 



P 





o 


+ c.cjg.p 


( 2 ) 


where A is the vector potential, e is the polarization vector, and q is the 

wave vector for incoming optical radiation, e is the magnitude of the 

charge of the electron, m^ is the free space electron mass, and p is the 

momentum vector of the electron in the crystal* The first term in (2) 

* abs 

gives the absorption, , and the second term gives the emission of 

photon. 

* 

Then, for a given interaction potential H , the transition rate from 

14 

the initial state ^ to the final state ^ for the absorption is given by 

W fi " I* I < *f I H lp bS |*i > I 2 « < E f “ E i - 

where E^ and E^ are the energies of the electron in the initial state and 
the final state, respectively and oj is the angular frequency of the inci- 
dent photon. If we neglect the interaction between the electrons in the 
well, the wave functions for the initial state and the final state 
after absorption can be written as^ 

*i “ u c (r) 5i (r) 

* A ^ 2 u c (r) e iK t’ r t (z), |z| < y 


(4a) 
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and 


* f * V < r > C f 

- A u c > (r) e i ^t* r t (z), |z| < -j (4b) 

-► + 

where A is the area of the well, L is the width of the well, k t , k t are the 

wave vectors of the electron in the x-y plane for the initial and the final 

states, respectively, r t is the position vector in the x-y plane, and 

u and u * are the cell periodic functions near the conduction band ex- 
c c 

tremum. The envelope functions ^ and satisfy the following 
Schrodinger f s equation in the effective mass approximation, ^ ^ 

2 2 

“ * ( z ) + l e l Fz $ < z ) * E <)> (z). |z| <. y » (5) 

dz 

and are given by the linear combination of two independent Airy functions 
Ai(n) and Bi(n), where r\ is defined by 

n ■ - t — 2 - m * v l 1 ^ 3 (E - |e| Fz) (6) 

(efiF; 

In Eqs. (5) and (6), m* and F denote the effective mass of an electron and 
the electric field, respectively. 

For intrasubband transitions, the matrix element |H op [ can 

be approximated by ^ 



5 


<♦, 


op 


»1> ■ <?f l«“l «i > 


eA .♦ ♦ <. 

2m~ q ’ r *'PUi> 
o 


2m 


e • <«f Ip I 5i> 


(7) 


21R 


CE i - E f ) t • <5, 


|r| C t >. 


where the portion of the cell periodic function has been taken care of 

following ref. 16, and we have used the dipole approximation. Since 

, L/2 

/ *i ^ * ^fi ' we that t ^ le ‘“J 01 ' contribution to the 

optical matrix element is the z-component of the r vector* So the absorp- 
tion is strongly polarization dependent. The absorption constant a in the 
well is defined, as^ 

, number of transitions per unit volume per unit time 
a ^ incident power per unit area 


rill l 

i f £ £, 

t t 


( 8 ) 


4dW 


fi 


i 2 a 2 \ 

K“» A o 


2yC 


where the summations over i and f are for the quantized initial and final 

energies respectively for the z-components of the momenta. If we calculate 

14 

the total transition rate and take into account the line broadening , we 
obtain 


! 


6 


- n 


y cm*k.gT e 


if 


2 2 
cos 0 |M fl | An 



(e - E (z) \ 

1+exp 

5 f s t 

V / 


L . E (^V 

1+exp 

F f 

i V / 


(r/2) 


(*„, - E fl ) 2 +(r/2) 2 


(9) 


with the matrix element 


M 


fi 


(z) (z) 

® (E, - E ) l/2 

— U/2 *V*> 2 *1 <2) d2 


( 10 ) 


(*) (z) Cz) (2) 

where E^ » E^ - E^ , and and E^ denote the quantized energy 

levels for the initial state and final state, respectively, p is the perme 
ability, c is the speed of light in free space, is the Boltzmann 
constant, T is the temperature, 0 is the angle between the polarization 
vector and the normal to the quantum well, n f is the refractive index, 

Rp is the Fermi energy which depends on the density of electrons in the 
well, and r is the line width* 

The oscillator strength f is given by 11 , 


2m <E< 2) - E (z) ) 2 

f - - i <z > 2 

•ir 


2 |m 


fi 1 


•fun^(Ej Z) - E< z) ) 


(ID 


In the zero field limit, f * 14*45 for the 1+2 transition, which is inde- 
pendent of the width of the well for the infinite potential well model* 

The experimental result*of f for L » 65 A (or * 101*27 A) is 12.2 

and slightly depends on the well width* 
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We calculated a for the first three states with field dependence 
numerically for T 3 300 K. In Fig. 2, we plot the absorption coefficient a 
for the incident photon with polarization perpendicular to the well (0 * 0) 
taking into account the first three states as a function of the energy of 
the photon with F ■ 0 (dashed line) and F * 250 kV/cm (solid line) for an 
effective well width 101.27 A, which gives the same ground state energy for 
F » 0 with the true well width L * 65 A and the barrier height a£ c » 245 

meV.** We use E * 6.49 meV which corresponds to about 2.5 x 10^/cm^ 

* 11 
electrons and r * 10 meV from the experimental results and is assumed to 

be independent for the variation of F. One can easily see that the tran- 
sition i 2 is dominant for both electric fields F a 0 and 250 kV/cm. For 
F * 250 kV/cm the absorption peak is shifted by 16 meV from 165 meV to 
181 meV and the peak amplitude is increased from 3153 cm * to 461 9 cm * • 

There are two distinct features for the case of the intrasubband absorption 

9 

compared with the exciton absorption. 

(i) The absorption peak for intrasubband optical absorption is 
increased in energy with increasing electric field over a wide range of the 
electric field, because for increasing electric fields the energy of the 
ground state decreases rapidly, while those of the higher subband states 
increase slightly then decrease slowly as cited in Ref. 6. On the other 
hand, for the exciton absorption, both the ground states of the electrons 
and the holes decreases. Thus the absorption peak is decreased in energy 
with increasing electric field. 

(ii) The absorption peak for intrasubband optical absorption i3 
increased in magnitude with increasing electric field because the electrons 



8 


are shifted to the same side of the well for both the initial and the final 

states with increasing electric field, and the energy difference E 2 - Ej 

also increases for the reason mentioned in (i)* As a result, the absolute 

value of the overlap integral for 1 + 2 transition increases* While 

for the exciton absorption, increasing electric field causes further 

separation of electrons and holes in the well and the decrease of the 

energy difference between the electron and the hole ground states* Thus, 

4 

the decrease of the absolute value of the overlap integral* 

It is also remarkable that the forbidden transition 1+3 for F ■ 0 
becomes possible when F is nonzero because the parity which prohibits the 
transition 1 + 3 no longer exists when F is nonzero. In Fig* 3 we plot 

as a function of F for the 1 + 2 transition, where is the 

value of for the 1 + 2 transition with F 31 0, One can easily see as 
expected that the ratio increases slightly from 1 as F increases. In our 
calculation, we assume r is constant, however, to account for the effect 
of the electric field on the absorption completely, further analysis of the 
electric field dependence of the line broadening is desired. 

In conclusion, we have calculated the electric field dependence of the 
intrasubband absorption within a conduction band of a quantum well. It is 
found that the absorption peak is shifted in energy and is also increased 
in magnitude with increasing electric field. The forbidden transition 
1+3 when F ■ 0 becomes allowable for the nonzero electric field. 

This work was partially supported by the NASA Grant No. NAG 1-500 and 

the Air Force Contract No. F33615-84-K-1557 . D. Ahn would also like to 
thank GTE for financial support. 
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Fig. 1 
Fig. 2 


Fig. 3 


Figure Captions 


. Potential energy profile for an infinite quantum well with width L 
subject to an external electric field F in the presence of 
incoming radiation with angular frequency /lioj. 

. Comparison of the intrasubband absorption coefficient a for infi- 
nite well with L ■ 101.27 A for the first three states with 
electron density 2.5 x lO^-'/cm^ electrons for the zero electric 
field (dashed line) and for the electric field of 250 kV/cm (solid 
line). 

. The normalized overlap integral fM 2 i/M2i ^ | 2 , where is for the 

zero electric field, is plotted versus electric field F. 
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: A Coupled Mode Formulation by Reciprocity and a 
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| Abstract— k eoapM mode formulation for parallel dielectric wave- 
uldes is presented via two methods: a reciprocity theorem and a vari- 
tional principle. In the first method, a generalised reciprocity relation 
|>r two sets of field solutions (£"’, ff (>> ) and (£ <J> , H m ) satisfying Max- 
ell’s equations and the boundary conditions in two different media 
' '(x, y) and y), respectively, is derived. Based on the generalised 

eciprocity theorem, we then formulate the coupled mode equations, 
j'he second method using a variational principle is also presented for 
general waveguide system which can be lossy. The results of the vari- 
jtional principle can also be shown to be identical to those from the 
eciprocity theorem. The exact relations governing the “conventional” 
|nd the new coupling coefficients are derived. It is shown analytically 
hat our formulation satisfies the reciprocity theorem and power con- 
ervation exactly, while the conventional theory violates the power con- 
ervation and reciprocity theorem by as much as 55 percent and the 
jlardy-Streifer theory by 0.033 percent, for example. 

I. Introduction 

T HE COUPLED mode theory has been very useful in 
the fields of integrated optics, semiconductor laser ar- 
rays or microstrip coupled transmission lines. A “con- 
ventional” coupled mode theory usually makes use of a 
perturbation theory to calculate the coupling coefficients 
[1], [2]. It has been recognized that a simple power con- 
servation argument for the powers in individual wave- 
guides leads to the fact that the two coupling coefficients 
and are complex conjugate of each other, which 
is generally not true if the guides are not identical [3]. A 
more rigorous approach has been recently proposed and 
very good numerical results have also been presented [3]- 
[3]. However, there is still considerable ambiguity about 
the reciprocity and the power conservation in the coupled 
mode theory. One knows that both the reciprocity relation 
and the power conservation are the two basic laws which 

E rnst be obeyed and they are usually used in electro mag- 
etics as necessary conditions to check the numerical ac- 
curacy [6], [7] of the results. The reciprocity relation is 
applied to the fileds and is applicable to a lossy medium. 
Thus, most results derived from the reciprocity relations 
do not contain any complex conjugate quantities. If the 
medium is lossless, the complex conjugate of the permit- 
tivity 6* equals to e itself, one then applies the conjugate 
fields to the reciprocity relation. On the other hand, the 
power conservation deals with the power and, thus, the 
complex conjugate quantities are usually used. 

Manuscript received April 11. 1986; revised June 4, 1986. This work 
wss partially supported by NASA gram NAG 1-500. 

_ Th e *«thor is with the Department of Electrical and Computer Engi- 
mering, university of Illinois at Urbana-Champaign, Urbana. 1L 61801. 
IEEE Log Number 861 1063. 


The goal of this paper is to present new coupled mode 
equations and analytical relations for the coupling coeffi- 
cients which follow the reciprocity theorem in a general 
lossy medium, and then the power conservation law if the 
medium becomes lossless. This new formulation removes 
the slight discrepancies of the power conservation en- 
countered in a previous theory presented in [3], [5]. The 
analytical relation governing the coupling coefficients Kob 
and is derived from a reciprocity relation for the fields 
instead of the power conservation law for the intensity. 
Thus, it is also applicable to any lossy (or gain) wave- 
guide system. 

The general reciprocity relation and the derivation of 
the new coupled mode equations are presented in Section 
II-A. A variational principle for a general lossy or lossless 
medium is presented in Section II-B while a previous 
method is limited to a lossless system [8]. We show that 
our formulation using the variational principle is identical 
to that of the formulation based on the reciprocity rela- 
tion. In Section HI, we derive the relation between the 
coupling coefficients and the propagation constants used 
in the coupled-mode equations. Note that this derivation 
is independent of the procedure in which one calculates 
those coupling coefficients and the propagation constants. 
We also show that the coupling coefficients and the prop- 
agation constants derived in Section II-A and Section 
II-B for the coupled mode equations do satisfy the reci- 
procity relation analytically. For a lossless case, the power 
conservation relation is derived from the reciprocity re- 
lation also. Finally, we present some numerical results 
and compare them with those of the previous theories. It 
is also demonstrated that an error of 55 percent in the 
power conservation using a previous theory [2] can occur 
unless the overlap integrals C M are taken into account 
properly. An error of 0.033 percent occurs using the 
Hardy-Streifer theory [3]-[5]. It is noted that the Hardy- 
Streifer theory, the theory of Haus et al. [8], and the 
present one give numerical results almost indistinguisha- 
ble on the plots of propagation constants and coupling 
coefficients for the examples considered so far, although 
slight differences exist among the three theories. 

II. Formulation 

A. Coupled Mode Theory from a Generalized 
Reciprocity Theorem 

In this section, we present a “generalized” reciprocity 
theorem for two sets of solutions (£ (1) , H iX) ) and (E a \ 
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H a) ) to Maxwell’s equations in nvo media e (1) and « (2> re- 
spectively. Based on the generalized reciprocity theorem, 
we show by choosing various e u> and e , and their cor- 
responding field solutions, two exact relations for the con- 
ventional coupling coefficients Kab and K ^ for two wave- 
guides a and b can be derived in Cases A1 and A2. We 
then derive the coupled-mode equations in Cases A3 and 
A4. A different approach using the variational principle 
for waveguide systems will be presented in Section II-B, 
and identical results of the two approaches are also illus- 
trated. 

1. A Generalized Reciprocity Theorem for Two Media 
f a> (x, y) and e m (x, y): Consider the first two Maxwell’s 
equations in a medium e w (x, y) 

V x £ (I) = 
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fxff 


( 1 ) 


— i«e ( ’£ 


(i)p(D 



(la) 

(lb) 


where the fields (E 0> , H (i) ) satisfy all the Maxwell’s equa- 
tions and the boundary conditions in the medium e 0) (x. 
y). For a different medium e (2) (x, y), the fields (£ <2) , H (2) ) 
satisfy a similar set of equations and also the boundary 
conditions in e (2) . Following similar procedures for the 
Lorentz reciprocity theorem, we obtain 

V • (£ ,(1) x - E m x H 0) ) 

= i«*® - e 0> )£ (1) • E a \ (2) 

If we apply the above relation to an infinitesimal section 
Az of a cylindrical geometry which is translational in- 
variant in the z direction, we obtain 

j j (£ (1) x H m - E {2) x ff (1 >) • tdxdy 

= i« | j (e a> (x, y) - e 0) (x, y))E (n • E a) dx dy 

(3) 

where the divergence theorem has been used. A similar 
equation using the polarization veqtor has been derived 
before [1]. However, our interpretation using e (2) and e (1) 
instead of the polarization vector is slightly different and 
will be shown to be very useful. We note that the above 
relations are exact as long as the fields (£°\ H (l) ) satisfy 
the Maxwell equations and all the boundary conditions in 
the medium e ll) (jc, y) and (£ <2) , H (1) ) in the medium i (2) (x, 
y), respectively. The above reciprocity relation is appli- 
cable to any two reciprocal media and is exact, while most 
reciprocity relations are applied to only one reciprocal 
medium with a polarization vector introduced and ap- 
proximated using a perturbation approach. The advantage 
of using the above exact relation will be shown in the next 
few cases when applied to a coupled- waveguide system. 
The time convention exp (-ioit) will be adopted in this 
paper. 

Case AJ: We choose first 




Fig. 1. .Schematic diagrams for various media under consideration: (a) 
y) with a single waveguide a . (b) y ) with a single waveguide 

b. (c) e(x, y) with both waveguides a and b. 


where e <a) (x, y) is a single waveguide a as shown in Fig. 
1(a), and we choose the solutions to be a guided mode 
propagating in the +z direction 

£<1> _ £<‘>+( X( y) = (£(«> + ££<«>) e ‘B.t (5a) 

H d) = n ia)+ (x, y) = (if< a) + mf) e i6ti . (5b) 

We then choose 


« 0) (x, y) * e ta, (x, y) 


(b) ( 


( 6 ) 


and 


£0) = £<*)- (x> y) e -'A« a (£<*> _ ££<*>) e ~^ 

(7a) 

H a) = H {b) -(x, y) = (-H\ b) + 2ff< w ) e~ i6 “ 

(7b) 

which are the guided modes propagating in the —z direc- 
tion for another waveguide b as shown in Fig. 1(b). Sub- 
stituting the above two sets of solutions into the reciproc- 
ity relation (3), we obtain 


« (l) (*. y) = e M (x, y) 


(4) 


Kto- K ak = i (Cab + Cba)0 b - 0a) 


( 8 ) 
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1 where 

' 

I K„ = ^ J j Ae^E, 0 ” • Ef - EfE^) dx dy 

“OB 

! (9) 

I » 

I = x ■ tdxdy (10) 

“OB 

I 

| and 

| Ae (,) (x, >') = e(jr, y) - e (,) (x, y), q * a, b 

(ID 

which are defined almost identically to those used in [3] 
except a constant factor of 4. The choice of the back- 
[ ground «(x, y ) is not unique (in general). Here, it is chc- 
| sen to be the coupled waveguide system (Fig. 1(c)) for 
convenience. One notes that K pq 's are the “conventional” 

I coupling coefficients except for the z components in the 
1 last term of the integrand [2]. Note that Equation (8) is 
an exact relation as long as the field solutions for each 
I waveguide system e <a> (x, y) and e (b) (x, y) are exact. For a 
j slab waveguide structure, the exact solutions are known 
and the identity (8) can also be proved analytically since 
| all the quantities K pq and C M can be derived. That proof 
is mathematically laborious but straightforward, and will 
not be shown here. Equation (8) is also a very useful re- 
flation in checking the numerical accuracies of the “ cou- 
pling coefficients” in_the computer program. One sees 
j clearly that in general when 0 b * 0„. Equation 

(8) shows the precise relation that the difference between 
the coupling coefficients is equal to the difference between 
the two propagation constants multiplied by the average 
of the overlap integrals and C*,. In the limit of ex- 
tremely weak coupling C&, « 1, we have AT* * 

Kgi, which is the reciprocity relation under the very weak 
coupling condition in a conventional analysis. 


Case A2: We choose 

« 0) (*, y) - e (a) (x, y) (12) 

E 0) =, E (a)+ (x, y) e iBtl (13a) 

H w m y) e ‘^ (13b) 

and 

e C) (x, y) - e {b \x, y) (14) 

E a) = E tb) +(x, y) e ,B “ (15a) 

fl® = fl< fc > + (x, y) (15b) 


Both fields are propagating in the +z direction. We obtain 
i from (3) 

I + KU - (K^ + K z *) = i (C ta - C*X/3* + 0 a ) 


where 


K> = " ! 

* 4 . 

| J Ae (,) E} p> • dx dy (17a) 

and 


(j) 

*--4. 

J J At^E^’E^ dx dy. (17b) 

Note from (9): 

K„ = K l M - K z „. (17c) 

(If only TE modes are excited, we have E\ p) = 0, p = a. 


b\ thus, = 0. Equations (8) and (16) wjll lead toCjafr, 
= CrifrQ In general, relations between and K'^, or 
Kit, and can also be derived from (8) and (16). In the 
following cases, we apply the reciprocity relation to the 
coupiea wavegude medium e(x, y) as shown in Fig. 1(c), 
and derive the new coupled mode equations. 

Core A3: We choose 

e (1) (*. y) = €(*, y) ( 18 ) 

and 

E™ * a(z) E? )+ (x, y) + Hz) E}‘ )+ (x, y) (19a) 

H™ = a(z ) H } 8)+ (x, y) + Hz) H\ b) + (x, y) (19b) 

for the transverse components. The above relations are 
just the modal expansions in terms of the two guided 
modes in waveguides a and b. We also note that the above 
expansion is only an approximate set of solutions to the 
Maxwell equations in the coupled-waveguide medium e(x, 
y) and the radiation mode has been neglected. Both wave- 
guides a and b are assumed to support only a single TE 
(or TM) mode. The extension to a multiple mode wave- 
guide is straightforward by including a summation over 
all the guided modes in each waveguide. The longitudinal 
components of the fields follow from Maxwell’s equa- 
tions for the waveguides 

.(«) ,<*) 

E< 0 a a(z) — E ( “ } (x, y) + Hz) — E[ b) (x, y) (20a) 

H[ l) » a(z)Hf(x, y) + Hz)H?'(x , y). (20b) 

A derivation of the above two components in (20a) and 
(20b) is given in Appendix A. A similar relation has been 
given for the z-component of the polarization vector in 
[11, and used in [3]-[5]. The factors e <a) l( and e (b) /e in 
(20a) have been ignored in [8]. We think they should be 
kept for consistency with the Maxwell equations as shown 
in Appendix A. 

For the second set of solutions, we choose the medium 
for a single waveguide a 

( a \x, y) = c^ 8) (x, y) 


(16) 


( 21 ) 
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and the guided mode solutions in the -z direction 

E a) = E (a) ~(x, y ) e~ iffmt (22a) 

H a) * H ia) ~{x, y) (22b) 


We obtain from (3) 

da(z) 4- Cjg db(z) 

dz 2 dz 


i(A + O a(z) 


* i (e. c ** c *'' + *») «a ( 23 ) 


where 

A, - j j J A« ( «W • £?» - ~ £^£?>) rfy. 


(24) 

To keep the same convention as in K„, the order p, q for 
the definition of A? is reversed from that in [3], and for 
later use. It is straightforward to show that A, satisfies 
the same relation (8) as by observing that 

A. = K„ + - A j J Ae ( «> EtfEP dx dy (25) 

where the second term is symmetrical when we exchange 
p and q. Thus. 


R„ - - K„ - K 9 » i (C„ + C»)(/3, - (3,) (26) 

which are exact relations. 

Case A4: We choose 


where the matrix elements for C and S are 


C » C * C v 

V v* '-'v -s 


where 


p,q = a,b (or 1, 2) 

(31) 

3* - *, + 8, (S'-fk) 

(32a) 

“ A» + Cpq&q 

(32b) 


and where and A* are defined in (10) and (24), re- 
spectively, and (26) has been used in (32). Wejiote that 
C, i * Cxi * C„ * Cj 2 ** 1, and the matrix C is sym- 
metric. The matrix S M is obviously symmetric following 
(32). 

Let 


— _ Q, + 
c — • 

We invert the matrix C and obtain the coupled mode 
equations 


da 

— = iy a a + iK^b 
db 

— - 


(33a) 

(33b) 


where 

7a - A + [A. + 03. - A)? 2 - £*c]/( 1 - c 2 ) 


« (1) (x, y) * e(x, y) 

and (£!'’, If) 0 ) and ff* l) ) to be the same as in the 
first set of solutions (19)— (20) in Case A3. We use for the 


second set of solutions 


« 0) (x, y) * e (b \x, y) 

(27) 

£® * E ib) ~(x, y) e~ i6kt 

(28a) 

H a) m H <b) ~{x, y) e~ i0kt . 

(28b) 

We obtain again from (3) 



- A + [A. ~ &.?]/(! ~ C 2 ) 

(34a) 

7* - A + [A* + (A - A)? 2 - ^c)/(l 

- c 2 ) 

- A + (As ~ - c 2 ) 

(34b) 

Kat = [A. + (A ” A - *»)?]/(! ~ C 2 ) 

* (A* - A**)/(l * C 2 ) 

(34c) 

A. - (As + (A - A - 0?]/U - ? 2 ) 

= (A. - A.?)/(i - c 2 ) 

(34d) 


± Q. da(z) db{z) 

2 dz dz 

- « (a — y ~ + As) a(z) + i(A + A*) b(z). 

(29) 

2. Coupled Mode Equations: Based on the results in 
Cases A3 and A4, we obtain the coupled mode equations 

®sH-*H ( 30 ) 

dz Lh(z)J Lb(z)J J 


where the first fonn in each equation is to compare with 
that in [3], and the second form is simplified after making 
use of (26) or (32b). One should know that although the 
matrices C and S are both symmetric, C~'S is not sym- 
metric in general. That is, Kai> * A*, unless we have two 
identical waveguides. This does not violate the reciproc- 
ity theorem or the power conservation law as will be pre- 
sented rigorously later. 

B. Coupled Mode Theory From a Variational Principle 
Applicable to a Lossy Medium 
Variational principle has been widely used to study the 
resonators, the waveguides or scattering from objects [9], 
[101. A general variational formula for the propagation 



! 
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constant 7 of the coupled waveguide system t(x, y) can 

to be (Appendix A): 


Je derived from two oppositely traveling modes of the 
System 

Jo) Jb ) 

£ 4 = af -j- £< a)± + at y E? H 

(42a) 

V, x £ + - = —iyt x £ + 

(35a) 

ff 4 = afH^ + a£H[ b)± 

(42b) 

V, x H * + = -i 7 2 x ff + 

(35b) 

where aj\ a£ are, in general, independent of af and d 2 . 
The variational formula can be put in a quotient of two 

V, X £~ — — iyi X £" 

(36a) 

quadratic forms 


V, X H~ + i(ixE~ ■ iyi x ff~. 

Dot multiplying (35a) by H~ and (35b) by E~ 
Ing the results, we obtain 

l /• 

(36b) 
, and add- 


(43) 





J J [E~ • (V, x Jf + + itaeE*) + W • (V, x £ + - iw^H*)) dx dy 


^ j j (H~ x £ + + £ - x ff + ) • 4 dx dy 


N 

D 


(37) 


where N and D denote, respectively, the numerator and 
he denominator in (37). A similar form to (37) has been 
ierived in [8], [10]-[12], except that we keep H~, £ + , 
stc., in the denominator which will be shown to be nec- 
essary in the following case. 

I It is straightforward to show that (37) is a variational 
formula for the propagation constant y by taking the first 
variation in y, 6y, from the trial fields 


where p, q * 1,2 or a, b in a two-waveguide system. 
The matrix elements C M are defined as (31). The deriva- 
tion of the matrix elements Q„'s is more complicated and 
is given by 

\ J J fwAe^- • Er 


£* 
h r* 


£o + SB* 


(38a) 

H$ + 8H 1 (38b) 

where Eq and H$ are assumed to be the exact solutions, 
(hat is, using (33) and (36) for ££ H $ and 70, one finds 


— k 


8 N 


A> J D 0 


[8N - y 0 8D] - 0 


(39) 


where N 0 and D 0 are the expressions in (37) evaluated 
tising £0 and H$. Thus any deviations of first order in 
}E ± and 8H ± only result in errors of second order (5£ 4 ) 2 
and (iff 4 ) 2 in 7. 

We choose the trial functions to be 

E; m a ,+£<»>* + <, + £<»>+ (40a) 

H; * affff + a 2 + ff<*> + (40b) 

[or the transverse components of the fields propagating in 
he +z direction. Here the subscript 1 refen to a, or 
waveguide a, and 2 for b for convenience. We also choose 

! E~ =* «,-£«- + of £{* >_ (41a) 

Hr * a;H\ a) ~ + afff<*>- (41b) 

for the transveree components of the fields propagating in 
the -z direction. The longitudinal components are found 


- • t x H? )+ 


- • i x £, W)+ 


+ iH^~ • V, x 



(44) 


where various relations such as those in Appendix A have 
been used. Using some vector identity and integration by 
parts for the last term in (44): 

j J fffW" ‘ * (* ““ £$»>+) dx dy 

- j j m w ^ Ef- £<« )+ dx dy. (45) 


We simplify Q M 

ft, - j J J • g?> - y EfE? 

• dxdy + y Cg ) g, 

= Rpq + Cpqdq 

» S„. (46) 


Thus it is clear that Q M is identical to S n . 
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The propagation constant y of the supermode is deter- 
mined from the variational formula (37). We thus take the 
partial derivative with respect to a~ regarding the ampli- 
tudes of the positive traveling waves a* to be independent 
of a~ 


and obtain 


frcs'.P _ 0 

dal 


S C„a* = ZQ„c 

<1 9 


(47) 

(48) 


where we have made use of (43) again. 
Noting that 


d_ 

dz 


hr 


for the system mode, we obtain the coupled-mode equa- 
tion 



-<> 2 - 



(49) 


which is identical to the coupled-mode equation (30) de- 
rived in Section Q-A since Q - S. If one takes partial 
derivative with respect to a* in (47), one obtains identical 
results as (49) since both C and Q are symmetric. 


ID. Reciprocity and Power Conservation 

Almost all the previous theories use the power conser- 
vation to find the relation between the two coupling coef- 
ficients K& and K^. That would lead to erroneous results 
since K& * in general, if two waveguides are not 
identical. An approximate theory from a more rigorous 
approach indicates some clue to the reciprocity relation 
by a power conservation argument but the results still 
contain some small discrepancies [3]. The explanation 
given in [3] was that they are due to the neglect of the 
radiation modes. In this section, we show that reciprocal 
relations can indeed be satisfied analytically and the pre- 
cise analytical relation of K+ to K* can be obtained, and 
the radiation field can be ignored from the beginning. The 
relation derived here should be obeyed and our coupled 
mode theory does satisfy this analytical relation. 

A. Reciprocity Relations 

Let us apply the reciprocity relation (2) to the two 
waveguide system described by t(x, y) (fig. 2). We 
choose « (1) (x, y) = e <2) (x, y) » e(x, y), and the two sets 
of solutions 

E ( , l) * a ( '\z) E? )+ (jc, y) + b m (z) £<* )+ (jc, y) (50a) 

B\" - * (1) U) *T<*. y) + b"\z) »?>+<*, y) (50b) 

and 

E? ) = c < 2) (2) £{ a) "(x, y) + b a> (z) E< b) -( x, y) (51a) 

= a®( z ) Hf~(x, y) + b m (z) H< b) -(x, y) (51b) 





Fig. 2. Two parallel dielectric waveguides applied to the reciprocity re- 
lation. The surfaces 5 t and Sj are normal to the £ -direction. The side 
surface Sj expands to infinity. The two sets of solutiosn used are: 1) 
« <n (-J) * 0, h m (-/) » K 0 , o o> (0) • equation (55a), b 0) ( 0) * equation 
(55b). 2) a <2) (-/) — equation (60a), » equation (60b), a (2) (G) 

- U 0 , b i2) { 0) - 0. 


for the transverse components where the radiation mode 
has been neglected. 

The volume of integration is chosen to be bound by 5,, 
Sj, and S d as shown in Fig. 2. Using the divergence theo- 
rem and the fact that the surface integral on the side S d 
goes to zero because of the radiation condition, we obtain 

j J (E\ l) x ff® - £p> X jr<») -tdxdy 

S\ 


- j j (E® X ff® - £® x ffj'*) - tdxdy (52) 

z-0 

Si 

which leads to 

a (1) ( 0) a 0) ( 0) + * Cba [<j®(0) b®(0) + a m (0) * a) (0)] 

+ h®(0) b a) ( 0) 


- a (,) (-/) a®(-l) + £± ±£s [„<•>(_/) b®(-/) 

+ a®(-i ) b w (-l )] + b (l \-l) b a \-l). 


(53) 

We next consider these two sets of solutions to be the 
coupled mode equations with two boundary conditions 
satisfied respectively. One starts at z » -/ with the 
boundary conditions 


(1) 


a®(-/) - 0 
b n) (-l) * V 0 


(54a) 

(54b) 



'HUANG: COUPLED MODE FORMULATION BY RECIPROCITY AND VARIATIONAL PRINCIPLE 
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I fl (,) (0) = V 0 i sin Me* 1 (55a) 


, n 2 ( x) = «<x)/« n 
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I-*, -I 


/,<»«)) = K 0 |cos M + i - sin (55b) 


7* - 7« 


| * = VA 2 + (57) 

(58 ) 

: The next set of solutions arc the propagating modes in 
the -i direction with the boundary conditions 

2) o®(0) = l/ 0 (59a) 

| b p) (0) * 0 (59b) 

ind the solutions when the mode propagates to z = — / 

a a \-l) = t/ 0 cos \pl - i ^ sin M e’* 1 (60a) 

I h (2) (-0 - C/ 0 — sin («0b) 

| y 

Substituting these field amplitudes (54), (55), (59), and 
(60) into the reciprocal relation (53), we obtain immedi- 
ately the relation 

! K ia -K ab = A (C* + C^,) (61) 


t A n, » A< lol («)/«, 






t— '3-1 


Fig. 3. (a) (b) (c) An illustrative example to show the two coupled wave- 
guides under consideration. There also exists an external perturbation 
between the two waveguides. 


plex conjugate fields are also solutions. We choose 
£,®(x, y) = a*(z) E\ a)1 (x, y) + b*(z) E (b)± (x, y) 

(63a) 

Efix, y) = a*\z) Hf'ix, y) + b*(z) H) b>1 (x, y) 

(63b) 


Vhich is the reciprocal relation that must be obeyed. Note 
fiat the above relation is exact and there is no complex 
ionjugate operation involved here. It is applicable to lossy 
js well as lossless systems. Each quantity in (61) can be 
Complex in general. Using our theory as derived in Sec- 
jton II, the quantities given by (34a)-(34d) do indeed sat- 
isfy the above reciprocal relation (61) analytically! The 
froof is straightforward by substitutions and making use 
if (26) for and Interestingly, the above relation 
61) is of the same form as (26) except that the propaga- 
ion constants are the modified y„ and y b instead of and 
!(, for individual waveguides. 

1. Power Conservations 
We choose the first set of solutions to be 
E?\x, y) = a(z) £, <fl)+ (x, y) + b(z) E* )+ (x, y) (62a) 

H?\x, y) = a(z) H\ a)+ (x, y) + b(z) H\ b)+ (x, y) (62b) 
or « (1) (x, y) = e (x, y). 

For the second set of solutions, we choose e (2) (x, y) = 
*(*» y). Since the medium is lossless, e* = t, the corn- 


making use of the z-inversion symmetry also. Substituting 
(62) and (63) into (52), we obtain 

P{z = -/) - P(z = 0) (64a) 

where 

P(z) = |a(z )| 2 + |b(z )| 2 

+ (Cab + QJ Re (a(z) b*{z)) (64b) 

turns out to be the power guided by the two waveguides, 
where we have used relations such as 


»J 

j £, (W x H\ a) • 

tdxdy 



j E (b) x Hf 

• i dx dy 

(65) 


for a lossless system assuming one chooses £, and H, to 
be real, which is possible [13]. We note that since the 
distance / between the two surfaces S ( and S 2 is arbitrary, 
(64a) leads to the fact that P(z) should be constant inde- 
pendent of z, which is also obvious from the power con- 
servation point of view. Using the boundary conditions 
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PROPAGATION 
CONSTANT (REAL PART) 
U/jmp) 



PROPAGATION 

CONSTANT (IMAGINARY PART) 
(l/pm) 



Fig. 4. The propagation constant* for the coupled waveguides in Fig. 3: 
(a) the real parts, and (b) the imaginary parts of the propagation constants 
(1/pm) are plotted versus the thickness (fim) of waveguide b. The exact 
solution (solid line), our resuits (dashed line), and the results using [3] 
(dotted line) are almost on top of each other. The crosses are results 
using [2]. 


that a(0) = 0 and b{ 0) = K 0 , we find [3] 


P(z) 


I K>| 2 |^1 + “jf [(Kab - Kba) 

+ ACCat + CfcJ] sin 2 tzl = 


Thus the “power conservation violation factor” for ex- 
citation in waveguide b at z * 0 

jp, ss 

* b-*a 


constant. (66) 


y l(Kab -K ba ) + MC* + C*)] (67) 
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COUPLING 

COEFFICIENTS 



<*) 


COUPLING 

COEFFICIENTS 



Fig. 5. (a) The real parts, and (b) the imaginary pans of the coupling coef- 
ficients and for the waveguide system in Fig. 3 are plotted versus 
the thickness (pm) of waveguide b. Our results (dashed lines) and the 
results using [3] (dotted lines) are on top of each other. The crosses are 
the results using [2]. 


phould be zero. One sees clearly that this condition has 
been derived in the previous section using the reciprocity 
theorem which is more general for lossy as well as loss- 
less cases. In deriving (67), one needs to restrict every 
quantity in (67) to be real for a lossless medium. Our new 
formulation presented in the previous sections does sat- 
isfy exactly these reciprocity conditions and power con- 
servation, since the factor F is zero if we substitute all 
quantities in (34a)-(34d) into (67) and use (26). The fac- 


tor F is an indication of the power conservation and the 
reciprocity relation. It can be used for the final numerical 
check of the consistency of the theory. Similarly, for an 
initial excitation in waveguide a at z = 0, one can define 
another factor 

Fa-t, = [<*. - **) - MC* + QJ] (68) 
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POWER CONSERVATION 
VIOL ATlONifjAC TOR (%) 
JO r ' ' « 


-,60 




* I 


05 r- 


H20 


F a F *••• * 1 

I u \^ ■ I • • ## ,*•••••••••••••••• • 


- - r. .-T r.-.--rn° 


H-20 


-.0Jh 


H-40 


.10 .12 .14 .16 .IS .20 

WAVEGUIDE THICKNESS l/xm) » 4 

Fig. 6. The power conservation violator factors f.-« and F.. k for the 
duee methods; our results (dashed line), the results using [3] (dotted 
lines) are shown using the left scale. The results using [2] (crosses) are 
shown using the right scale. 


to check the numerical accuracy. The numerical results of 
these two factors in (67) and (68) using various methods 
will be presented in the next section. 

IV. Numerical Results and Discussions 

In this section, we consider an example from [3]. The 
coupled mode equations (33a) and (33b) with the expres- 
sions in (34a)-(34d) are used in the numerical calcula- 
tions in this paper. The refractive index profile is shown 
in Fig. 3 where an external perturbation between the two 
planar waveguides also exists. We choose the index vari- 
ation to be along the x direction, and for TE polarized 
waves, the electric field has only the y-component. The 
refractive indices are n t = n 5 = Re (n } ) = 3.4, n 2 = n 4 
= 3.6 and an additional loss exists between two guides 
such that n| - n] = il. 299 x 10" 3 . The other parameters 
are t 2 * 0. 15 fim, r 3 = 0.4 fim, X = 0.8 fim, and r 4 varies 
from 0. 1 fim to 0.2 fim. The numerical results using an 
exact root-searching approach have also been shown as 
die solid lines in Fig. 4(a) and 4(b) for the real and imag- 
inary parts of the propagation constants. (In [3], the “ex- 
act” numerical method combines a root searching ap- 
proach assuming a lossless system to find the real parts of 
the propagation constants, and a perturbational approach 
for the imaginary parts when the loss is added. The final 
“exact” results in [3] are indeed very good compared with 
our exact root searching approach.) The results of the the- 
ory in this papeT are shown as the dashed lines, and the 
results using that in [3] are shown as the dotted lines. We 
see clearly that all three methods agree very well with 
each other. The results of a conventional method [2] are 
also shown as the crosses which deviate more from the 
exact solutions especially for the imaginary parts of the 
propagation constants. 


In Fig. 5(a) and 5(b) we compare both the real and the 
imaginary parts of the coupling coefficients using our 
method and the methods in [2] and [3]. It is clear that our 
results do agree very well with those using the method in 
[3] with a different approach, which has been checked with 
the “exact” numerical results presented in [3]. We note 
that our results satisfy the reciprocity and power conser- 
vation analytically and, thus, the factors F b -. a and F a -i 
in (67) and (68) are zero while the F s of the method in 
[3] still contain a small discrepancy which is around 0.033 
percent at a maximum value at r 4 = 0.1 fim, and the F’s 
of the method in [2] yield a maximum power discrepancy 
of 55 percent at r 4 = 0.1 fim (instead of only about 20 
percent as claimed in [3]). Detailed calculations of the 
two power conservation violation factors are shown in 
Fig. 6 (assuming the lossless case, i.e., Im Inf - nf} = 
0) where F fc _ fl and for the method in [3] are shown 
(die dotted lines) in the left scale. The results are within 
0.033 percent. The results using [2] (crosses) show in the 
right scale that F a - k for excitation in waveguide a has an 
error of power conservation of 21 percent at r 4 = 0.1 fim, 
and F b ~ a for excitation in waveguide b has a value of 55 
percent. Our results (the dashed line) for F b -. a and F a ^ b 
are always zero or within the round off errors in the com- 
puter, and the power conservation is indeed satisfied. 

One should also note that the relations using the reci- 
procity and power conservation laws are necessary con- 
ditions, not sufficient conditions, for the accuracy of the 
numerical results [7]. They usually serve as checks, not 
direct proofs, of the numerical solutions to the Maxwell 
equations. 

V. Conclusions 

A new coupled mode formulation has been described 
via two methods: a generalized reciprocity relation and a 
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l variational principle. Both give the same results. Exact 
| analytical relations governing the coupling coefficients 
i Rat, and (also and K^) and the propagation con- 
stants of individual waveguide (3 a and /3» (y a and y b ) are 
i derived. These relations are used to show that our for- 
mulation does satisfy the reciprocity theorem and the 
power conservation analytically. Numerical results com- 
pared with the exact solutions and a previous method [3] 
which contains a slight discrepancy show that our new 
formulation should be very useful and self-consistent. We 
hope this paper will also clarify the reciprocity relation 
for the coupled waveguides. 

Appendix A 

Derivation of (20a) and (20b) 

For the guided modes, we have 

V, x E (fl)+ - iunH w * « X E w+ (Al) 

i 

I V, x H (a) * + iue w E ia) * - -i& a i x H (e) * (A2) 

and a similar set of equations for t (b) , £ (t)+ , H (b)+ , and 
i For the coupled-waveguide medium, we have 

| V, X E — iufiH * —iyi x E (A3) 

| V, X H + iutE * -iyi X H. (A4) 

Breaking the equation into the transverse and longitudinal 
I components, we have 

I E t - -4— V,xH, = -4- (a(2)V, x Hl e) + 

— JO* 

+ b(z)V, x Hl b) +) 

Jo) Jb) 

= 0 ( 2 ) — Ef* + biz) — E (b) * (AS) 

which is (20a) ih the text. A similar procedure can be 
applied to H t and leads to (20b). 
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A Coupled-Mode Theory for Multiwaveguide 
Systems Satisfying the Reciprocity Theorem 
and Power Conservation 

SHUN-UEN CHUANG 


Abstract — Two sets of coupled-mode equations for multiwaveguide 
systems are derived using a generalized reciprocity relation; one set 
for a lossless system and the other for a general lossy or lossless system. 
The second set of equations also reduces to those of the first set in the 
lossless case under the condition that the transverse field components 
are chosen to be real. 

Analytical relations between the coupling coefficients are shown and 
applied to the coupling of mode equations. It is shown analytically that 
our results satisfy exactly both the reciprocity theorem and power con- 
servation. New orthogonal relations between the supermodes are de- 
rived in matrix form with the overlap integrals taken into account. 

I. Introduction 

T HE COUPLING of mode theory in parallel wave- 
guide systems has been of great interest in applica- 
tions to directional couplers, laser arrays, waveguide 
switches, etc. [1], [2]. Although it has long been recog- 
nized that the previous coupled-mode theory is only ap- 
plicable to very weakly coupled systems I3]-[8], signifi- 
cant improvements for strongly coupled waveguides have 
only been presented recently in series of papers [5], [8]- 
[ 11 ]. 

The major improvement is probably the inclusion of the 
overlap integrals C pq defined in [8], when evaluating the 
power, and its resultant corrections to the various param- 
eters such as the propagation constants and the coupling 
coefficients in the coupled-mode equations. Using two dif- 
ferent methods, one based on a generalized reciprocity 
theorem and the other based on the variational principle, 
a new set of coupled-mode equations has been derived for 
a general (lossy or lossless) system [12]. Both methods 
give the same results. 

In this paper, we apply the generalized reciprocity theo- 
rem [12] to a multiwaveguide system. The lossless case 
is treated here separately from the general lossy case, since 
in a lossless system, one may prefer to deal directly with 
powers for which the complex conjugates of the fields are 
needed, while for the general lossy case, one may not re- 
quire any complex conjugate operations in the formula- 
tion [8]-[12]. Thus, the definitions for the overlap inte- 
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grals and the coupling coefficients presented in Section III 
will be different from those for the general lossy case pre- 
sented in Section IV. As will be shown in this paper, only 
when one chooses the transverse electric and magnetic 
field components to be real functions, the two formula- 
tions will be identical in the lossless limit. New properties 
of our coupled-mode equations are also presented analyt- 
ically with the overlap integrals properly included. 

n. Generalized Reciprocity Relation 

Assuming that the electric and the magnetic fields E ( ", 
H (,) satisfy the Maxwell equations in a medium e < 1 > ( jc, 
y) (for the whole space) and the corresponding boundary 
conditions and that E <2> and H <2) satisfy the Maxwell 
equations in another medium e <2> (x, y) and the corre- 
sponding boundary conditions, it is straightforward to 
show that [12], [13] 

' V • (£<'> x H <2) - £ <2> x H u> ) 

= / w (t <2> -e <,> )£ <,> • £ <:> (1) 

with the same procedure used for deriving the Lorentz 
reciprocity relation [14], When applied to a cylindrical 
geometry with an infinitesimal distance in the z-direction, 

(1) reduces to 

^ j j (£<’> X H <2) - E <2> x H ny ) ■ tdxdy 

= Ua j j (€ <2 > (z, y ) - e <0 (z, y))E <l> • E <2> dx dy . 

( 2 ) 

Here e ( 1 2 (x, y ) and e ( 2 > (or, y) can be general media such 
as a single waveguide or a multiple waveguide system as 
long as they are translational invariant in the z-direction. 
The time convention exp ( -iut) will be used in this pa- 
per. One notes that the two reciprocal relations (1) and 

(2) are exact as long as the two sets of field expressions 
(£ <l> , H a> ) and (E <2> , H <2> ) are exact solutions to 
the Maxwell equations in medium 1 > (x, y) and e ( 2 > (x, 
y) respectively. 

HI. Coupled-Mode Theory for a Lossless 
Multiwaveguide System 

In this section, we derive the coupled-mode equations 
for a lossless multiwaveguide system. 
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A. General Properties of the Fields of the Guided 
Modes 

When the medium e(x, y) is lossless and translational 
invariant in the z -direction, one knows that the field so- 
lutions of the form exist 

a) (E, + E z ) e iBz 
(H, + H z ) e iBz 

which correspond to the fields propagating in the +z di- 
rection. Here, we assume the above set of solutions to be 
the guided mode of the system. Based on inversion sym- 
metry in the -z-direction, the following set of fields will 
also be solutions to the Maxwell equations [8], [13]— [15] 

b) (E, - E z ) e' iBz 
(-H, + H z )e' iBz 

which correspond to the fields propagating in the -z-di- 
rection. If the medium is lossless e*(x, y) = e(x, y ), by 
taking the complex conjugate of the Maxwell equations 
or applying the time-reversal concept, it is easy to show 
that the following two sets of solutions also exist: 

c) (E? + Et) e~ iB ' z 

( -H * - H*) e' iB ‘ z 

d) (E* - £*) e iB ' z 
W - H?) e‘ B "' 

where the * sign means complex conjugate. Since we con- 
sider the guided modes of the lossless system (excluding 
the leaky modes, cutoff modes, etc.), the propagation 
constant 0* is real. It is thus clear from a) and d) that one 
can choose the transverse field E, to be real, and find im- 
mediately that H, is real; E z and H : are purely imaginary. 
However, if one uses complex E, (e.g. , in an optical fiber 
with a circular cross section, £,( p, <t>) can be of the form 
J m (k p p) «"*•), one finds that H„ E z and H z will also be 
complex. From these general properties of the field so- 
lutions, we next derive the coupled-mode equations for a 
lossless system and some analytical relations between the 
coupling coefficients and the overlap integrals. 

B. The Derivation of the Coupled-Mode Equations 
CASE (1): Suppose we choose 


H 

II 

(3) 

£<>> = (£<«> + £<«>) e iB " z 

(4a) 

£ n> = ( H\ q) + H[ q) )e^ z 

(4b) 


to be the guided mode propagating in the + z-direction in 
a medium e (,) (jc, y) with a single waveguide q. We also 
choose for the second set 

e <2> (x, y) = e ip) '(x, y) = y) (5) 


and 

£< 2 > - (E^* + £<*)*) e~ iB ’ z (6a) 

H <2> = (- H \ p) ’ - H[ p) ’) e~ iB ' z (6b) 

which are also solutions as discussed before. They cor- 
respond to the fields propagating in the — z-direction. 
Substituting the two sets of expressions into the general- 
ized reciprocity relation (2), we obtain 

- r‘„ - (g, - g,) ( c ” j c * ) (7) | 

where 

oo 

K„ m ^ J J &e iq) (E] p)m • E (q) + E[ p), E (q) ) dx dy (8) 

| 

and 

• i 

C„ = ^ J J E (q) x H\ p) * • tdxdy. (9) ! 

— 00 

We note that (7) is an exact relation since the fields (£* 1 \ i 
£ 0) ) and (E <2> , H <2) ) are exact solutions to Max- 
well’s equations in « <i> (jc, >’) and e <2> (x, y), respec- 
tively. 1 

CASE (2): In this case, we choose e< 1 J (x, y ) to be the 
medium of the multiwaveguide system e(x, y) 

< <1> (x, y) = «(x, y). (10) 

The solutions to the system are given approximately by 

£,<■> * S a ? (z)£|”(x,y) (11a) 

p*\ 

£ ( <0 = Yj a q {z)H\ q \x,y) (lib) 

p _ i 

for the transverse field components. The z-components are 
given by 

£<» = S^^(:)£{'»(x,y) (11c) 

N 

£<» * S a,(z)£<«>(x,y). (lid) 

p m i 

A similar derivation for the above relations has been given 
in [13] for the polarization vector or in [12]. One notes 
that E^Cx, y), q - 1, • • • , N are not orthogonal func- 
tions, and the overlap integrals C p<J * 0. The second set 
of solutions is chosen as 

€ <2> (x, y) = ( ip) ‘(x, y) = e lp> (x, y) (12) 

E< 2> « (Ej pr + Ei'>*) e~*' z (13a) 

ff< 2 > = (-Hj pr - H[ p) ‘) e~ iBpz . (13b) 

Substituting the two sets of expressions (10)-(13) into the 


I 


t 

I 


I 

I 


I 
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generalized reciprocity relation (2), we obtain 

2£„|a,(z)-iE(** +/3/„K(z) (14) 
where 

0D 

k M = J J J Af ( ’^£, (p) * • £{*> + — £< p) *£<«^ dx dy 

(15) 

and 

C*^(C„ + C*). (16) 

One notes that £ p? is a hermitian matrix, £ M = £ J,. It is 
straightforward to show that satisfies the same relation 
(7) as Kpq because 

k„ = + j J j Ae <4) A« (p) £j P) *£j ,) dr dy (17) 

where the second term is equal to its complex conjugate 
quantity if one exchanges p and q where both Ac (p> and 
A f (?) are real (lossless). Thus 

J? — k* = j? _ j?* 

-</>,- ■(«,-«, >e„ 

(18) 

which is an exact relation. It is seen clearly that only if 
8 P = £ ? , one has (or if the overlap integrals 

are very small in the extremely weak coupling case, k M 
k^). Otherwise one should treat k M and K*, as differ- 
ent quantities in general. One defines the matrix elements: 

4* - + js /„ - 4 + £„$,. ( 19) 

Thus, the coupling of mode equations can be written as 


- * 

(20) 

where £ is clearly hermitian since 


“ &/>9 

(21) 

which can be shown from (19) and a is a vector with its 
elements given by a q (z), q = 1,2, • • • , N. Another 
way to write the above equation is either 

e ^ - .(<* + *<?)« 

dz 

(22) 

or 


£ -j- = i(k + dB)a 
dz 

(23) 


where J? is a diagonal matrix with the elements given by 


the propagation constants of individual waveguides & p . 
Here the superscript + means complex conjugate and 
transpose of the matrix. The second form (23) is useful 
since 


— = i(B + C~'/t)a (24) 

dz 

while the first form (22), which is similar (but not iden- 
tical) to that of [10], [1 1], requires more algebraic manip- 
ulations in evaluating (£"'b£ + d~ l £* ). 


C. Power Conservation 

In Section III-B, we derived the coupled-mode equa- 
tions in matrix form (20), where £ is related to the overlap 
integrals C M and C£, and £ is defined in (19). Both £ 
and Q are proved to be hermitian without any approxi- 
mation in the matrix elements. Let us look at the power 
guided along the muitiwaveguide system 

P(z) - | Re j J E, x H* • t dx dy 
* Re I Ea*(z) C M a,(z) 

L P'4 

= Z a*(z) a f (z) = a + (z) da(z) 

(25) 

where £ is defined in (16). If the medium is lossless, the 
power of the guided mode must be independent of the 
position z, i.e., dP /dz — 0. We have the lossless con- 
dition 





da^(z) 

dz 


0. 

(26) 


Using the coupled-mode equation (20), one finds imme- 
diately that the lossless condition is equivalent to 

/£<(&, -$*)«,« 0 (27) 

where we have used the fact that £ M is hermitian. Since 
a* and a, can be arbitrary values, we obtain 

= 0. (28) 

That is, must be hermitian, which is true, since we 
have shown that is indeed hermitian in (21) from the 
definition (19). Thus our formulation satisfies exactly the 
power conservation. An example will be shown later 
which illustrates this power conservation criterion. 


D. Power Orthogonality of the Supermodes 

Let us choose two sets of solutions to be two distinct 
supermodes in the multiwaveguide system e(x, y) 

e a> (x, y) - e(x, y) (29) 
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(30a) 

//, <0 = ^ 2 a^H^ e h,z (30b) 

where a {,) with elements a^\ q = 1,2, * * • , N, is the 
eigenvector for the first supermode with a propagation 
constant 7 , 

e <2> (*, y) = e(x, y) = «*(*, y) (31) 

E< 2 > = (^2 a^'E\ p) '^j e-** (32a) 

ff< 2 > = ^-^2 e~ lyji (32b) 

and a ij) with elements a ( p J \ p = 1 , 2 , * • * , N, is the 
eigenvector for the second supermode with a propagation 
constant 7 j. The reciprocity relation (2) gives 

( 7 , - 7; ) 22 a^ m {C M + C*) = 0. (33) 

Since 7 , * 7 , we obtain the general orthogonality condi- 
tion: 

d 1)+ £a U) = 0, j. (34) 

That is, any two eigenvectors corresponding to different 
propagation constants are orthogonal to each other with a 
weighting matrix given by t. 

An alternative way of deriving (34) is simply by look- 
ing at the coupled-mode equation (20). The supermode 
solution a(z) is given by the form 

a(z) = ae iyz . (35) 


Thus, the matrix equation (20) for the coupled-mode 
equations reduces to the eigenequation 

7 to = fa. (36) 

The eigenvalue y satisfies 

det | (5 - 7 C | = 0. (37) 

Since both <? and (5 are hermitian, the eigenvalues for (36) 
must be real, that can be shown from elementary matrix 
theory [16]. It is also obviously true from the fact that the 
medium is lossless. Another property of the matrix equa- 
tion (36) is that two distinct eigenvectors a u) and a (J) are 
orthogonal to each other with the “weighting matrix” d 

a U) *6a {i) = 0 , i * j. (38) 

One notes that in the extremely weak coupling case, the 
coupling of mode equations have the same form as ( 20 ) 
except that £ should be replaced by /, the identity matrix. 
Thus, the orthogonality relation (38) reduces to the well- 
known results: a u> *a u> = 0 for i * j in conventional 
theory. 


IV. Coupled-Mode Theory for a General (Lossy 
or Lossless) Multiwaveguide System 
In general, a multiwaveguide system can be lossy. The 
previous formulation will not be applicable anymore. Ac- 
tually, the formulation for a lossy medium is very similar 
to the formulation in the previous section, except one does 
not have any complex conjugate operation, and special 
care is taken for the ^-components of the fields as can be 
seen from Section IE- A. A derivation has been presented 
in [12] which is similar to that in the previous section. 
Therefore, we briefly give the results below. 

A. The Derivation of the Coupled Mode Equations for a 
Lossy System 

CASE (1): Following the procedure in Case (1) of Sec- 
tion ffl-B, except that we choose the second set of solu- 
tions to be the form b) in Section III-A, it is easy to derive 
[ 12 ] 

K„-K qp = (f3 p - 0 ,) — j — ( 39 ) 

where 

OD 

K„ = ^ j J A« ( «>(E! p> • £ ^ - E[ p) E { f') dx dy (40) 

— 00 

and 

00 

c * = H \ £<,,> x h<,p) ■ f ^ (4i) 


where no complex conjugate operation is involved, and 
there is a negative sign in the integrand of (40). The above 
definitions (40) and (41) are the same as tho$e used in [ 8 ] 
except for the constant factor of 4. The difference is only 
apparent because once we choose the normalization con- 
dition C,| * C 22 - • ‘ • C NN - 1 , the factor of 4 is ab- 
sorbed in E, and H,. Thus, numerically, K M is identical 
to that in [ 8 ]. 

CASE (2): Following the procedures in Case (2) of 
Section m-B, we choose the first medium and the field 
solutions to be the same as ( 10 ) and ( 11 ), and the second 
medium and the field solutions to be 


e <2 >(*,y) = < <p) (*,y) (42) 

£ <2 > * (El p) - £< p) ) e~ i9 ' z (43a) 
ff< 2 > = {-H\ p) + H[ p) ) e- iB " : . (43b) 

We obtain again from the generalized reciprocity relation 

. 2 ) 






+ « C ” + C " 

<? ' 


0„(z) 


2 


(44) 
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where 

K„ = ~ j j Ae'^E^ • £<*> - ~ £<'>£<«>) (45) 

and C pq has been defined in (41). One can also show that 
K pq satisfies the same equation as K pq in (39) by recogniz- 
ing that 

K„ - ~ J J Ae<’> £<'>£<*> dx dy (46) 


where the second term is symmetric with respect to p and 
q. Thus 



K pq -K qp = ( s p - e q ) ~— 2 c - 

(47) 

which is also an exact relation. 

The coupied-mode equation can be written in a 
form 

matrix 


c^«(z) = iQa 

(48) 

where 

T + C„ 

2 

(49) 

is symmetric, and 



Cpq + C^ 

Q„ = k 9 + 6, - - 

c„ + c v 

~~ Kpq + 2 ~ QdP 

(50) 

is also symmetric. The matrix equation can also be writ- 
ten as 


d 

— a(z) = iMa 

(51) 


m = r'g 

(52) 

or 

M = C~'BC + C~'K t 

(53) 

or 

M = B + C~'K 

(54) 


where B is again a diagonal matrix with the propagation 
constants 0 P as the elements, and the superscript T means 
transpose of the matrix. Equation (53) is compared with 
the form in [ 10],_[1 1] . One sees that the only difference 
is that the matrix C is used here while the matrix C is used 
in [8]-[lI]. AT r used in this paper is the same as £in [10], 
[11] from the definition (40) except for the factor of 4. 
The final form (54) is simpler than (53) since B is used 
instead of C~ l BC. Thus, our coupled-mode equation 
looks simpler using the form (51) with M given by (54), 
than that in [10], [11]. The coupling coefficients K pq de- 
fined in (15) for the lossless case or K M defined in (45) 


for the general case differ from K„ q (defined in (8)) or K pq 
(defined in (40)) by the factor e (p ’/( in the second part of 
the integrand. This factor is also taken as one in [17]. We 


believe that it should be more self-consistent to keep the 
factor since it was derived from (11c) making use of Max- 
well’s equation as shown in [12, appendix A]. 

B. General Orthogonality Property of the Supermodes 

Following a similar procedure to that in Section III-D, 
one applies the reciprocity relation (2) to any two super- 
modes 


£, <0 = ( 2 e iyiZ 

(55a) 


»<'> ■ a^H^ e iy,z 

(55b) 

and 

£, <2> = ( 2 d p j) E\ p) ^J e~ iy>z 

(56a) 


tf, <2> = (- 2 a ( p i) H\ p) \ e-‘ yiZ . 

(56b) 

Using the eigenvector a U) = column ( a\'\ a 2 \ 
a^) and a similar form for a u> , one obtains 



a <J>r Ca U) = 0, for y, * y } 

(57) 


which is the reciprocity relation that should be satisfied 
by any two eigenvectors of the matrix equation 

yCa = Qa (58) 

which follows (48). Alternatively, because both C and Q 
are symmetric matrices, the general orthogonality relation 
(57) is a well-known property in matrix theory [16]. 

C. Reciprocity Relation for Two Sets of Solutions with 
Separate Boundary Conditions 

Let us look at the boundary value problem for a set of 
solutions to the coupled-mode equation (48). The general 
solution a(z) is given by 

a(z) = Ae ,Tz A~ l a(0) (59) 

where 

T = diagonal (y,. y 2 , • • • , y N ) (60) 

for a given boundary condition a(0) and the wave prop- 
agating in the +z direction. Here the matrix A is defined 
to have the ith column given by the ith eigenvector of the 
matrix (58), and >,,•••, y N are the eigenvalues of (58). 

Consider a first set of solutions at z = 0 given the con- 
dition a a> (z = -/), 

«<'>(()) = Ae ,Tl A-'a n >(-/). (61) 

Let us look at another set of solutions with the boundary 
condition given at z = 0 and the wave propagating in the 
-z direction to z — —l 

a <2> (—l ) = Ae ,T 'A -l a <2> ( 0). 


(62) 
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Fig, 1. (a) A single waveguide p described by t ip) (x, y) in whole space, 
(b) A single waveguide q described by y) in whole space, (c) A 

multiwaveguide system described by e ( jt. y) in whole space. 


Applying the reciprocity relation (1) to a cylindrical sur- 
face enclosing the planes z - —l and z ~ 0 with a radius 
going to infinity, one finds [12] 


j j (£ u> x H w - £ <2> x fl <n ) ■ tdxdy 

Z = -/ 


= j j (£' <1> x H< 2> - E a> x H <1> ) • tdxdy. 

z~0 


(63) 

From the previous two sets of solutions, we have (Fig. 2) 
£, <u - 2 a< ,> (z)E ( t p> (64a) 

P m l 

= 2 a { p '\z)H\ p) (64b) 

p - 1 


and 


E{ 2) - 2 aq 2> (z)E ( , 9) (65a) 

9*1 

H, <2> - - 2 u' 2> (z)ff<’ ) (65b) 

9-1 

w here the second set of fields propagates in the -z-direc- 
l ion. The vector a (l> (z = 0) is related to a ( 1 > (z « -/) 
(61), and a <2> (z * — / ) is related to a <2> (z = 0) by 


o 


o 


I 


o; 

o' 



I 

I 


q 


N-l 


N 


;o 

:o 



i 


0 


2**1 


2*0 


Fig. 2. A multi waveguide system with possible excitation either at wave- 
guide p at z 88 —l and the wave propagating in the direction, or at 
waveguide q at z « 0 and the wave propagating in the -z-direction. 


(62). The reciprocity relation (62) reduces to 

« <1>T (— /) <7a <2> ( — /) = a<’ >r (0) Ca <2> (0) (66) 

or equivalently 

fl <1>T (-f) CAe ,T U-'a <2> (0) 

= o <1>r (— /) (Ae ,T 'A-') r Cfl <2> ( 0). (67) 

Since the initial conditions for a> ‘ ‘ ( -/ ) and a ( 2 > can be 
'0 


«<’>(-/) = 


»< 2 > 


( 0 )- 


0 

0 

0 

0 

1 

19J 


pth position 


(68) 


(69) 


-* qth position 


where p and q can oe arbitrarily set between 1 and N, we 
find the reciprocity condition: 

CAe ,TI A~' = (CAe^'A-') 7 (70) 

where we have used the fact that C T - C. Since the matrix 
A has each column given by the eigenvector of (58), we 
have 


CAT = QA 


(71) 


or 


r — A~ l C~ l QA = A~ l MA. (72) 

Substituting the above relation into (70), we find that the 
reciprocity condition (70) is the same as 

CM * (CM) 7 (reciprocity condition) (73) 



I 
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i.e., the product CM must be symmetric. We see clearly 
that our formulation (48) satisfies this condition because 

CM - Q (74) 

where Q has been proved to be symmetric in Section IV- 
A. To illustrate our results, we show in the next section 
the special cases for two coupled waveguides and three 
coupled waveguides. 

V. Special Cases 

In summary, the coupled-mode equation is put in ma- 
trix form 

c|a = iQa (75) 

where both C and Q are symmetric. In another form, it is 
given by 

j 

7 «* iMa (76) 

dz 

where 

M - C~'Q. (77) 

The reciprocity condition (74) requires CM (or Q) to be 
symmetric. The above formulation is very general and is 
applicable to both lossy as well as lossless systems. 


A. Two-Coupled Waveguides 
If N = 2, one has 



M = 

7a ^ab 



(78) 



~^ba 7 b - 



where 






7a 

= 3, + (*„ 

- C 12 *2,)/(l ~ 

C\ 2 ) 

(79a) 

7 b 

= 02 + (*22 

- C 12 * 12 )/(l - 

Cj 2 ) 

(79b) 


= (*i 2 ~ *22 

C,2)/(l ~ 

• C] 2 ) 


(79c) 

Kta 

= (*2I ~ *,, 

C 12 )/(l - 

' C 12 ). 


(79d) 

As has been ^pointed out in [8], the overlap 

1 integrals C 12 


and C 21 or C I2 are obtained from the integration over 
whole space in the transverse direction, and can be sig- 
nificantly large even if is small. Thus the factor 1 — 
C ] 2 may become very_small and *„* is large. The reci- 
procity condition that CM be symmetric gives 

K ab -K ba = (y a -y fc )C 12 (80) 


which has been shown in [12], and can also be proved by 
substituting (79a)-(79d) into (80). The two eigenvalues 
7 i, 7:. and eigenvectors are well known: 


I 

| 

7i - 

(81a) 

i 

72 = $ ~ ^ 

(81b) 

! where 

1 

, 7 b + 7a 

*■ 2 

(81c) 


^ = 

VA 2 + K+Kta 

(81d) 

A = 

7i - 7a 
2 

(81e) 

a <«> = 

[ K* 1 
La + +\ 

(82a) 

a (2) = 

[“«; *] 

(82b) 


where the orthogonality relation a a)T Ca (2) = 0 is indeed 
satisfied and it is the same as the reciprocity condition 
(80). 

B. Three-Coupled Waveguides 

UN = 3, we have 

M = B + C~'K (83) 

which can be calculated easily by inversion of C, noting 
that C is symmetric: 

C \2 — C 2 1 = (C l2 + C 2 j)/ 2 (84a) 

c, 3 = c 31 = (C, 3 + Cji)/2 (84b) 

C„ — Cji - C 33 = 1. (84c) 

The reciprocity condition that CM is symmetric leads to 
e.g .,(CM) n = (CM) 2l 

m \2 ~ m 21 ~ C\2 ( m 12 — m 22) + Ql m 31 — Q3 m 32 

(85) 

which will be useful later. Let us consider a symmetric 


case with the two outer waveguides identical [10]: 

K n = *32 * *21 = *23 ( 86a) 

*,3 = *3, (86b) 

*„ = *33 * * 22 - (86c) 

The matrix elements of M are obtained from (83) 

m ll = m 33 = 3l + [*|1 ( 1 — C j 2 ) — * 2 ,C 12 a — C 13 ) 

+ * 13 (C ? 2 - C n ))/D (87a) 

"*72 = ft + [* 22(1 + C. 3 ) ~ 2* 12 C 12 ] (1 - C n )/D 

(87b) 

m n ~ m 32 - (Kn ~ ^ 22 ^ 12 ) (1 — Cu)/D (87c) 
m 13 = m 3 l = [*13(1 — C) 2 ) — *2,C 12 (1 — C 13 ) 

+ *u(c? 2 - C 13 )]/D (87d) 

m 2 1 = m 22 = [^2l(l + Q 3 ) — (^11 

+ *,3)C,2](1 - C l3 )/D (87e) 

D = (1 - C I3 )(1 + C 13 - 2C] 2 ). <*lf) 
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The total guided power is given by | 

P(z) = Re [a + (z) Ca(z)] | 

= a + (z) Ca(z) = 1 + Fsin 2 y (94) , 
where the factor F is given by 

F = [m 12 (l + Cn) ~ m 2l 

+ C i2 (m22 - m u - mu)]. (95) ' 

For a lossless system, the power conservation requires that j 
P(z) be independent of the position z. Thus the factor F 
provides a check of the energy conservation. Although it 
has been shown before that our formulation satisfies the 
energy conservation exactly by using the fact that the two 
matrices C and Q are hermitian, one can also see from the | 
reciprocity relation (85) substituted into (95) that F is in- 
deed zero provided that we choose E\ p) and H\ p) to be ' 
real functions for a lossless system. Therefore, C pq 
Cpq = real. Numerical results will be given in the next 
section. The factor F from two previous methods [4], [10] 
will also be calculated. j 

VI. A Numerical Example and Discussions 

In this section, we illustrate our theory by a numerical 
example [10] and compare it with those of two previous 
methods [4], [10]. We consider three coupled waveguides 
with the two outside waveguides identical and symmetric j 
with respect to the center waveguide (Fig. 3). Using the j 
theoretical results discussed in Section V-B, we calculate j 
K m , C m , and & p . The analytical relation (47) is used to 
check the numerical accuracies of these quantities. We , 
show in Fig. 4(a) the three eigenvalues yj, y 2 , and y 3 j 
from (88a)-(88c), which are the propagation constants of j 
the three supermodes, versus the separation t between the 
waveguides. We compare our results (dotted line) and the 
exact solutions (solid lines) of the multilayered structure 
in Fig. 4(a) and those of the method in [10] (dashed lines), 
the method in [4] (crosses) in Fig. 4(b). We see clearly 
that the results using the method in [10] and our theory 
agree very well with the exact calculation. There is a slight 
error for the third eigenvalue y 3 near cutoff where the sep- 
aration t is reduced to near 0.2 /tin. In our calculation, we 
choose the same parameters as in [10], n = 3.4, n, = 

3.6, n 2 = 3.63, </, = d 2 = 0.15 fim , and t varies. The 
method of [4] clearly has larger errors in y, and y 3 , es- 
pecially y 3 deviates from the exact results over a wide 
range of t near cutoff. The result of y 2 using three methods 
agree with each other very well because 

72 = m u ~ m ii = + (X)i — K]s ) / (1 — C 13 ) 

(96) 

using (87) and (88). Since C| 3 and (AT n - K l2 ) are very 
small, (C, 3 = 0.136 - 0.00436 and (*,, - ff, 3 ) = 
-0.0237 0.0004 (1 / M m) at t = 0.2 /xm ~ 0.6 ^m). 
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(a) 






(C) 

Fig. 3. Three coupled waveguides under investigation, (a) e (x) / e 0 for the 
three-waveguide system, (b) At ll, (x)/e 0 . (c) 4f <2) U)/«o- 


if we have C , 3 = 0 (theory of [4]) the difference is neg- 
ligible. Since C 13 = C 31 = Ci 3 , the theory of { 10 ] gives 
the same results for 7 ; as the results of this paper. 

In Fig. 5, we show the power conservation violation 
factor F for the initial excitation at the center waveguide, 
«(0) = column (0, 1 , 0). We see clearly that our results 
indeed satisfy the power conservation very well and F is 
always zero. The factor F calculated from [10] is always 
very small (less than 0.08 percent), but F calculated from 
[4] can be as much as 42 percent. 

Numerically speaking, our results are as good as or 
slightly better than those obtained from [10], The new 
features are that our formulation is derived using a sim- 
pler approach, and it satisfies both the reciprocity theorem 
and the power conservation law analytically, while [ 8 ]— 
[ 11 ] can only show numerically that their method satisfies 
the power conservation and the reciprocity theorem ap- 
proximately. (One should note that power conservation 
and reciprocity are only satisfied self-consistently and not 
exactly since the modal expansions ( 11 ) of the fields are 
approximate.) Our formulation also leads to the general 
orthogonality relations (38) and (57) with the overlap in- 
tegrals properly taken into account, that cannot be ob- 
tained from the formulation in [ 8 ]-[llJ. By setting the 
matrix C or C to be the identity matrix, the coupled-mode 
equations and the orthogonality relations all reduce to the 
results of a conventional analysis [4], Our numerical re- 
sults also show that ignoring the overlap integrals does 




Fig. 4. (a) A comparison of the propagation constants of the three super- 
modes between the exact results (solid lines) and the results of this paper 
(dotted lines) for the three coupled waveguides in Fig. 3. ..(b) A compar- 
ison with other two methods: results using [10] (dashed lines), results 
using [4] (crosses). The results of this paper are given by dotted lines. 



Pig. 5. The power conservation violation factor F for an excitation at the 
center waveguide of Fig. 3. The results of F using (10) (dashed line), 
and results of this paper (dotted line) are shown using the left scale, and 
the results of [4] (crosses) follow the right scale. 
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lead to erroneous results violating the power conservation 
significantly as has also been pointed out in [8] and [12]. 

VII. Conclusions 

Two sets of coupled-mode equations for a multiwave- 
guide system have been derived using a generalized re- 
ciprocity relation, one set for a lossless, and the other for 
a lossy or lossless system. New general orthogonality re- 
lations between the eigenvectors of the supermodes have 
been derived. We have derived the conditions on the ma- 
trix elements for the reciprocity theorem and the power 
conservation laws and have shown that our formulations 
do indeed satisfy those conditions analytically. 
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APPLICATION OF THE STRONGLY COUPLED-MODE THEORY 


TO INTEGRATED OPTICAL DEVICES 
S. L. Chuang 

University of Illinois at Urbana-Champaign 
Department of Electrical and Computer Engineering 
Urbana, Illinois 61801 

ABSTRACT 

A theory for strongly coupled waveguides is discussed and applied to two- 
and three-waveguide couplers and optical wavelength filters. This theory makes 
use of an exact analytical relation governing the coupling coefficients and the 
overlap integrals. It removes almost all of the constraints imposed by a 
simpler and approximate coupled-mode theory by Marcatili. It also satisfies the 
energy conservation and the reciprocity theorem self-consistently. We show very 
good numerical results with the overlap integral as large as 49 percent. The 
applications to electrooptical modulators, power dividers, power transfer devi- 
ces, and optical filters are all presented with numerical results. 
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1. INTRODUCTION 

The applications of the coupled-mode theory in integrated optical devices, 
such as waveguide couplers [1] - [3], laser arrays [4], [5], and optical filters 
[6] - [8], have been well known. However, theoretical improvements for strongly 
coupled waveguides have only been attempted very recently [9] - [14]. A simple 
and approximate version of coupled-mode equations for parallel dielectric wave- 
guides has also been presented by Marcatili [15] to account for the asymmetric 
properties of waveguides using a newly found relation between the coupling 
coefficients and the overlap integral of two coupled waveguides. A few con- 
ditions are assumed in that paper: 

(1) A scalar formulation of the fields is considered. 

(2) The refractive index perturbation is very small such that second-order 
terms can be ignored. 

n 2 * n 2 [l + A a (x,y) + ^(x^)] 2 

2 (l) 

* Uq[ 1 + 2A a (x,y) + 2A b (x,y)] 

Thus the new relation between the two coupling coefficients in [15] is only 
approximate. 

(3) The overlap integral c is assumed to be small (weakly coupling) and is 
not included in the coupled-mode equations because the coupled-mode 
equations in [15] are almost the same as those for the conventional 
theory [2] without including the overlap integrals in the four 
coupling parameters, y^, y^ 9 and K^ a . 

In this paper, we apply the theory developed in [9] - [14] and show that 
all the above conditions are not required. It is shown that an exact analytical 
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relation governing the coupling coefficients, the overlap integrals and the pro- 
pagation constants derived in [13] using a generalized reciprocity theorem can 
be combined with the formulation of Marcatili and will give very good numerical 
results even for strongly coupled waveguides. It has been pointed out in [9] 

that the four parameters Y » Y, > K and K should include the overlap 

a b ab ba 

integrals to obtain correct propagation constants of the supermodes. Since only 

^a* ^b individual waveguides and of the conventional coupling coef- 

ficients are used in the coupled-mode equations in [15], that theory will not 
yield accurate numerical results and may violate energy conservation signifi- 
cantly [9], [13] unless the overlap integral c « 1 which is assumed in [15]. 

In Section 2, we briefly review the strongly coupled-mode equations derived 
in [9] - [11], [13] and [14], their orthogonality relation, and an exact rela- 
tion between the coupling coefficients a ^d K ba # We a ^ so s ^ ow that this 

exact relation can also be derived from power conservation or reciprocity rela- 
tion for a lossless medium. In Section 3, we consider the two coupled wave- 
guides combining the strong coupling of mode equations and the formulation of 
[15] and illustrate the electrooptic effect. We then study the three waveguide 
couplers as power transfer devices and power dividers. Previous studies in [16] 
and [17] assume all three guides have the same refractive indices, and a direct 
numerical approach for the multilayered structure is taken. In Section 4, the 
cross-talk problems for both two and three coupled waveguides are investigated. 
Numerical results are presented and compared with those in [18] and [19]. In 
Section 5, the application of the coupled-mode theory to the optical wavelength 
filters is studied and the theoretical results are compared with the experimen- 
tal results in [8]. Finally, we give conclusions in Section 6. 
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2. THEORY OF STRONGLY COUPLED WAVEGUIDES 
Three very similar formulations of strongly coupled waveguides have been 
presented in [9], [12], [13]* The formulation by Haus et al. [12] is limited 
to the lossless system and has a small difference in the z-component of the 
electric field for the trial functions in the variational approach. The formu- 
lation of Hardy-Streif er [9] does not satisfy energy conservation and the 
reciprocity theorem and still contains a small error, while the theory of [13] 
(which was derived in a much simpler way) satisfies these laws analytically. 
Independently, a reformulation of [9] has been made [20] recently and is iden- 
tical to that of [13] after the modifications. The coupled-mode equation in 
[13], [14] and its properties are summarized below. 


2 . 1 Strongly Coupled-mode Equation 

The coupled-mode equation in vector form is given [13], [14], [20] by 


c — a(z) - Qa(z) 


( 2 ) 


or 


dz 


a(z) - Ma(z) 


(3) 


where 


Q - CB + K 

M - C -1 q - B + C _1 K 


(4) 

(5) 


where the vector a(z) has each element a^(z) given by the electric field ampli- 
tude for the transverse component of the mode in waveguide p, and the matrix 
elements C and K are defined in Appendix A. The matrix B is a diagonal 

pq pq rr 8 
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matrix with the diagonal elements given by g^, •••» 6^ of each individual 

waveguide in the absence of all other waveguides. It should be noted that (1) 
the two matrices C and Q are symmetric [13], [14], [20] and that is very important 
to prove the orthogonality property of the superraodes, and (2) the matrix M is 
not necessarily symmetric in general. 

2.2 Orthogonality of the Supermodes 

The supermodes of the multiwaveguide system satisfy the orthogonality 
relation [14] for symmetric matrices C and Q in Eq. (2) 


a (i)T Ca ( ^ - 0 


for * Yj 


( 6 ) 


where (a^^) is the eigenvector of the superraode with the propagation 

constant (y^), anc * su P erscr ^P c T denotes the transpose of the matrix or 
vector. 

2.3 Reciprocity Condition 

To satisfy the reciprocity relation, one finds that [14] 


CM - (CM) T 


(7) 


i.e., the matrix Q must be symmetric which is true as derived in [14]. For two 
coupled waveguides , 


M 


Y a K ab 
K ba Y b 


Equations (7) and (8) give 


C ab ' ^a * (Y a ‘ V< 


where c 3 C 


C 12 + C 21 


12 


( 8 ) 


( 9 ) 


2 


(see Appendix A). The above formulation is true in 
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general for both lossy and lossless systems. If the system is lossless, one may 
also have a slightly different formulation as presented in [12], [14] using 
field quantities involving complex conjugates. 

2 .4 Power Conservation 

If the multiwaveguide system is lossless, one can choose the transverse 

field components E and H to be real functions and find that E and H are 
t t z z 

purely imaginary and C and K are real [14]. The total power guided by the 
multiwaveguide system is 


P(z) »^-Re // E t x • zdxdy 


a (z)Ca(z) 


( 10 ) 


where the superscript + denotes the conjugate and transpose of the vector a(z), 

and one has chosen E^^ and to be real. Thus using the fact that C and Q 

dP ( z ) T 

are real matrices, one finds that the condition — — » 0 also leads to Q * Q 
— — T 

or CM * (CM) , which is the same as the reciprocity condition in (7). Actually, 
condition (7) is very general since it is applicable to both lossy and lossless 
cases. A similar formulation (for a lossless medium) leads to the fact that Q 

* m 

is Hermitian provided one uses complex conjugate quantities with and 

matrices as defined in [14]. The Hermitian matrix becomes obviously symmetric 

when it is real. Another derivation of the lossless condition for two coupled 

waveguides is shown in Appendix B, which also leads to Eq. (9) when y , y. , K , 

a b ab 

K, are real, 
ba 
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3. TWO AND THREE COUPLED WAVEGUIDES AND IMPROVEMENT OF 
MARCATILI f S THEORY 

In this section, we present a combination of the vector formulation for 
strongly coupled waveguides and Marcatili's theory which assumes two weakly 
coupled waveguides. We also discuss the applications to three coupled wave- 
guides used as either power transfer devices transferring power from one outer 
guide to another or as power dividers. The electrooptic effects when these 
devices are used as modulators are discussed. 

3.1 Two Coupled Waveguides 

3.1.1 Improvement of Marcatili's formulation 

We start with the coupled-mode equations 

If ’ l V“ + iK ab» <“*> 

£ • ‘V * 1K ba b <Ub > 


where 


Y a “ 8 1 + (K 11 " c 12 K 21 )/d _ c i 2 > 

(12a) 

— —2 


Yb * ^2 T ^22 ~ ^ 1 2^12 ) / (* " ^ 12 ^ 

(12b) 

K ab " <K 12 ‘ K 22 C 12 )/<1 ' C 12 ) 

(12c) 

K b. ' (!t 21 - *U=12 )/<1 ' 

(12d) 


where the subscript 1 refers to waveguide a or 1 , and 2 refers to waveguide b or 
2, whichever is convenient. 
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One notes that in the theory of Marcatili [15], (1) is assumed to be 
zero in the above four parameters, (2) and are ignored, and (3) 
and are defined only for scalar fields (pure TE case). Thus that formula- 
tion is almost the same as that for the conventional theory [2] and will lead to 
significant errors if becomes larger than, say, 10 percent (C^ and C ^ are 
normalized to be 1) [9], [13]. One notes that an exact relation holds between 
the conventional coupling coefficients [13], 


K 12 " K 21 


(s x - b 2 ) 


C 12 + C 21 


(13) 


while a similar relation found in [15] is only approximate since Che derivations 
there have assumed the refractive index variation & a (x,y) and A^(x,y) « 1 
(which is a good practical approximation). Using this relation, one can show 
that the following relation is true using (12), (13) and c - * ^1’ 


K ab ' \a ‘ l \ " V c 


(14) 


which is precisely the reciprocity condition, and it is the same as the power 

conservation condition for a lossless case (Appendix B). We define the asyn- 

chronism factor [15] in terms of the more correct parameters Y > Yi_» K ^ and 

a b ab 

K in (12a)-(12d). 


5 


Y b ~ Y a 

2/K K 
ab ba 


( 15 ) 


Given the initial excitation at z • 0 of a two-coupled waveguide, a(0) ■ 1, 
b(0) » 0, we obtain [9], [13] 



(16b) 


iK 

bU) - — — sin 

’I' 




where 


Y b + Y a 
♦ 2— 


(17a) 


# 

A 


y4 2 + K ,K, 

“ ab ba 

Y b ' Y a 
• 


(17b) 

(17c) 


It is easy to show also 





(18) 


or 



e 


sinh 



(19) 


The solutions (16a) and (16b) can be written as 


ail) - |cos A ab K ba Jl(l +« 2 ) l/2 


- i 


an 2 ) 1/2 


sin 


W (1+ « > 


2x1/21 




(20a) 


sinh _1 C(5 


b(£ ) - i 


sin [/K 7 k Tjl (1+5 ) ' ] e 


2vl/2, i$£ 


(20b) 
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The output power in waveguide a when waveguide b terminates at z - l is 
obtained using 

E t (x,y f z » l) * a( 2)E^ a ^(x,y) + b(£)E^ b \x,y) (21a) 

- I u^ a) E^ a)n (x,y) (21b) 

n»l 

H t (x,y,z = 2.) a a( 2)H^ a \x,y) + b( l)H^ b ^ (x,y) (22a) 

- I v^ a) H^ a)n (x,y) (22b) 

n**i n 


where the expansion in (21a) or (22a) is in terms of individual waveguide modes 

and in (21b) or (22b) is in terms of all the guided and radiation modes of wave- 

( a ) 

guide a alone since they form a complete set [9]* Multiplying (21) by and 

integrating over the cross section, one obtains 

u| a ^ =* a(Z) + • (23) 


Similarly, one finds 

v<*> - a(t) + C 21 b(i) . 


(24) 


These boundary conditions at z * 0 and z * l follow very closely those in 
[15]. The guided power due to the first mode 8^ in waveguide a is, thus, 


_ 1 „ r (a) (a)* 1 ff (a) 1 _(a) 1 . , 

p a ■ j Re ^ u i v i 2 i j E t x H t * zdxdy3 


, , 1 ■ C 12 C 21 \ 2sinh* l c6 , 2 r , ,2 x l/2, 

1 - ^ e sin l K. 2( 1 + 5 ) ] 


1 + a 


ab ba 


( 25 ) 
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using Eqs. (20), (23) and (24). A similar procedure for the output power in 
waveguide b when waveguide a is terminated at z = l leads to 


P fe = Re[(C 21 a + b)(C* 2 a* + b*)] 


C 12 C 21 + 


1 ' Cl2C21 sin 2 mi - *v /2 ) 


1 + 6 


(26) 


These results are very similar to those in [15] except the parameters are 
defined in terms of the more accurate parameters y , y, , K and K, . 

3 D 3. u u 3 

3.1.2 Numerical results for two strongly coupled waveguides 

In Figs. 1(a) - (d), we show numerical results for two coupled Ti-dif fused 

LiNbO^ channel waveguides modeled as two coupled slab waveguides (which is 

possible using the effective index method [11]) with the refractive index in 

waveguide a, n * 2.2 + , the effective refractive index in waveguide b, 

a t 

n, * 2.2 - -■—> where the refractive index difference 
d z 


An = n - n (27) 

a b 

is proportional to the externally applied voltage V across the two waveguides. 

The refractive index outside the two waveguides is assumed to be constant, 

n n = 2.19. The waveguide dimensions are d « d, =* 2 ym; the edge-to-edge 
u a b 

separation t * 1.9 ym. The wavelength X is 1.06 ym. In Fig. 1(a), the 
asynchronism 5 is plotted versus the refractive index difference. We see 
clearly that | 6 | is linearly proportional to | An } • The overlap Integrals 
C ^2 (dashed line) and (dotted line) with their arithmetic average c (solid 
line) are shown in Fig. 1(b), where they vary between 0.168 at An =* 0 to around 
0.178 which do not satisfy the condition in [15] for weak coupling (c < 0.1). 
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The numerical results for the propagation constants are calculated with more 
than 7 digits of accuracy, and the energy conservation law [14] is also checked 
to be valid with errors always less than 10 (In [14], two energy conser- 
vation violation factors have been defined and are used to check the numerical 
accuracy of the results,) Figure 1(b) differs slightly from the qualitative 
drawing of [15] with a range of variation in the overlap integrals around (0,178 

- 0.168)/0.168 » 6.0 percent. The coupling coefficients K , K, , and 

ab ba 

/K _K_ are shown in Fig. 1(c), which agree well with the qualitative results of 
ab ba 

[15], The output powers P (solid curve) and P. (dashed curve) are shown in 

a b 

Fig. 1(c). They do agree very well with the qualitative drawing of [15]. One 

notes that the minimum of P does not occur right at An = 0 (where P, * 1.0) due 

a b 

to the cross-talk problems which are discussed in Section 4- The power P^ 

actually goes to almost zero (P fl = 0.00051 - 33 dB) at An 2 -0.0002, where 

P^ reduces to 0,9723. The asymmetry of and the symmetric properties of 

P, versus An or the applied voltage agree very well with what has been presented 
b 

in [15] . However, our numerical approach provides very good numerical results 
even for the strongly coupled case with c > 0.1, while the theory of [15], 
although taking into account the asymmetry properties of coupled waveguides, is 
limited to weak coupling cases. 

3.2 Three Coupled Waveguides 

Let us consider a symmetric case for which the two outer waveguides are 
identical. Solutions for this case have been obtained in [10], [11], and [14] 
and will not be derived here. 

3.2.1 Power transfer devices 

When used as power transfer devices, the three coupled waveguides are 
assumed to have an initial excitation at z 3 0 



12 


a(0) 


and the input power is easily found to be 1. 

The solutions at position z are found to be [10], [11] 






sin U e 2 + e 1 ^ 


2 m. 


a 2 (z) 


21 


sin (|2) e 1 ^ 


and 


a^Cz) 


cos f - i —if 
2 ;|r 


sm 




i* 


2 -e^2 z 


where 


i|f a j "" ®22 ^ ^^12^21 

♦ * m H + m i3 + m 22 


and the three propagation constants of the supermodes are 


_ * + 
Y 1 T 


(28) 


(29a) 


(29b) 


(29c) 


(30a) 

(30b) 


(31a) 


y 2 * " m 13 


(31b) 



y 3 


(31c) 
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where the matrix elements have been derived in [10], [11], and [14]. 

The output power in waveguide 1 at z » £, where waveguides 2 and 3 termi- 
nate, is 


P out,l ■ * C l2 « 2 (t> ♦ C l3 a 3 (t)![3 1 Ct) * C 21 . 2 (l> 

following a similar procedure as in Eqs. (21) - (26). The output power at wave- 
guide 3 when waveguides 1 and 2 terminate at z * i is 


P 

out 


,3 


Re{ [C 31 a 1 C A) + c 32 a 2 (£) + + c 23 a 2 ()l) 

+ 3^ (O] } • 


(33) 


When applying (32) and (33) to a power transfer device, one may need to assume 
|Cj^ I and l C 3 jJ are stna11 since waveguides 1 and 3 are not terminated. 

3.2.2 Power dividers 

When used as power dividers, the three coupled waveguides have an initial 
excitation at z 3 0, 


a(0) 


0 

1 

0 


and the input power P can be found to be 1 . 


UO], [14] 


(34) 

The solutions at position z are 


2m 


^(z) 


a 3 ( Z ) - i- ^ 


12 


sin 






(35a) 


a 2 (z) 





(35b) 
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One finds the output powers in waveguide 1 and 3 to be equal using (35) in (32) 
or (33) since the two outer guides are identical. 

3.2.3 Numerical results for power transfer devices and power dividers 
The three coupled waveguides considered here are assumed to be symmetric 
with respect to the center guide. We assume the dimensions of the three wave- 
guides to be d^ d^ 38 d^ ■ 2 pm. The edge-to-edge separation of two nearby 
waveguides t * 1.9 pm. The wavelength \ is 1.06 pm. The refractive indices are 
assumed to be 

Oj - n 3 - 2.2 + (36a) 

n 2 = 2.2 - (36b) 

where the refractive index difference between either one of the outer guides and 
the center guide An =* n^ - n^ is proportional to the applied voltage V. We 
first plot the unnormalized asynchronism 2y^ - versus the refractive 

index difference An. One sees clearly at An a 0 the fact that all three guides 
are identical does not imply the synchronism condition 

2y 2 “ Tfj " Y 3 - 0 (37) 

is met. At An = 0, we find 

Y x - 13.0172261 
Y 2 * 13.0138696 
Y 3 = 13.0094738 

and 



0.0010393. 
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Choosing the coupling length l to be fixed at * 2 n/ where 

^0 * ^ Y 1 ~ y 3^ at An * 0 (38) 

we find the output powers P t (solid curve) and P _ 0 (dashed curve) as 
r r out,l out ,3 

shown in Fig, 2(b) when the waveguides are used as power transfer devices. Peak 
power transfer from guide 1 to guide 3 actually does not occur at An * 0 as can 
be seen from Fig, 2(b), This is due to the cross-talk problem when the syn- 
chronism condition is not met. It occurs actually at 2y 0 “ Y, “ Y^ = i.e., 

when An 2 -0,00023, The cross talks are calculated in Section 4. When used 
as power dividers, the three coupled waveguides are assumed to have a coupling 
length l - tt/iJJq * L^/2. The output powers in guides 1 and 3 versus the 
refractive index difference An are shown in Figure 2(c), Maximum output power 
does occur at An * 0 for the power dividers. 

In Figs, 3(a) and 3(b), we show the output powers P . and P - versus 

out , 1 out , J 

the coupling distance l normalized to - 2 tt/(y^ “ Y^) for An * 0 (L^ * L ^ 

in this case). Since the synchronism condition is not met, p out ^ does not go 

to zero due to the cross-talk problems. Both P . and P do not show 

periodic behaviors for the power transfer devices that have been discussed in 

[10], [11], although only the magnitudes of |a^(i)|, \ or |a^(jt)| instead 

of powers are presented there. This nonperiodic behavior is due to the 

asynchronism (2 y 2 “ Yj ~ Y^ * 0) of the three supermodes. When this condition 

is met (it occurs at An * -0.00023), P . and P _ 0 do show periodic beha- 

out,l out, 3 

viors as shown in Fig. 4(a), where * 2 tt/(Y| “ Yj) Is evaluated at that 
An 4= 0. However, the output powers for the power dividers always show periodic 
functions because of the symmetry of the excitation and the waveguide structure 
(only two of the three supermodes, Yj. anc * Y^> are excited). It is easy to see 


16 


from Eq. (35) that the output powers will be periodic functions of the distance 
l for both Figs. 3(b) (2y 2 - Yl - Y 3 * and (2y 2 ” Y 1 ~ Y 3 = °^* 
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4. CROSS-TALK PROBLEMS 

Cross-talk problems have been investigated recently for two coupled wave- 
guides [18] and three coupled waveguides [19] using either the conventional 
coupled-mode theory or direct numerical approach for the propagation constants. 
We apply the strongly coupled mode equations here to investigate the cross-talk 
problems* 

4*1 Cross Talks in Two Coupled Waveguides 

In the design of two coupled waveguides, one usually chooses the coupling 
length j i such that 

ipl m nn/ 2 , n - odd integer (39) 

and b(fc) is maximum, a (i) * 0 provided that A * 0, i.e., two waveguides are 
identical. One finds immediately that the output power in waveguide b is maxi- 
mum. However, the output power in waveguide a is not zero because there is 
still an overlap of fields between modes in waveguides a and b. This cross-talk 
power is easily obtained by setting 5 * 0 in Eq. (25) as a conservative estima- 
tion [18] 

— ? 

Extinction ratio » P^(6 m 0) » G i2 C 2 1 * C 12 (40) 

since C^ * ^21 w ^ en two waveguides are identical. The formula (40) only pro- 
vides a conservative estimation since it assumes waveguide a continues. In 
reality, guides a and b may start to separate at z * i gradually. Thus, (40) is 
only an approximation [18]. 

This result showing that the cross talk is proportional to the square of 
the overlap integral agrees with that obtained in [18]. However, our numerical 
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calculations show that for two-coupled GaAs waveguides with the dimensions 

d * d, 3 2 yin, the edge-to-edge separation t * 1.9 ym, the refractive indices 
a d 

n a n, ■ 3.44, and the outside refractive index n a * 3.436, the cross talk is 
ah u 

-10 dB, which is close to -12.6 dB of [18] but not identical. We believe our 
number here is more accurate since we have calculated the propagation constants 
and B 2 u P to 7 digits (after the decimal point) of accuracy; the power con- 
servation and the exact analytical relations are all checked so that the errors 
are always less than 10 The studies of cross talks in [18], [19] assume 
identical waveguides and the refractive indices are fixed (An » 0). One finds 
using Fig. 1(d) that the cross talk P can actually be given by Eq. (25). At 
An = -0.0002, the extinction ratio goes to zero! Thus a very good extinction 
ratio can be obtained with a slight asymmetry introduced in the two waveguides 
with An * 0 . 

4.2 Cross Talks in Three Coupled Waveguides 

Three coupled waveguides have been introduced to decrease the cross talks 
when used as power transfer devices from one outer waveguide to another. 

However, the synchronism condition 

2 Y2 " Yi ~ Y3 “ 0 

needs to be satisfied; otherwise, the cross talks may be proportional to the 
overlap integrals C^ and C 23 of the two nearby modes instead of the two outer 
guided modes C 13‘ When used as power transfer devices, one chooses 

l * 2tt / ( y ^ “ Y 3 ) “ 2 */^ (41) 

or an integral multiple of i such that a ^ (jO * 0 * We find the extinction ratio 


to be (using sin * 0 ) 


I 
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Extinction ratio « P 


out , 1 

2 


sin 


sin 


2 y 2 - Yi - y 3 

2 2 

4 - r 

fa - n * n " 

4 * 

H J 

' cos 

[_ 4 

> 2 ■ n " y 3 i 

. r 2 2 

'" 2 Y 2 - Yi - Y 3 ~] 

L 2( Yl • Y 3 > 'J 

+ cos 

- 2 <n - y 3 > "J 


(42) 


and similarly, the output power 

~ 2 y 2 ~ n - y 3 


P „ 0 * cos 
out,i 


2( *1 - Y 3 > 


2 2 

+ C 3 sin 


2 Y 2 ~ Yi - Y 3 
- 2(y 1 " *3 5 


(43) 


i This analytical result for cross talks is useful since it explains clearly 

(I) If the synchronism condition is met, 

| 

1 2 Y 2 ■ Yi - Y 3 - 0 


I 


Extinction ratio 


C 


2 

13 


j which is expected. 

(2) If the synchronism condition is not met, the first term will contri- 
j bute, and it will be proportional to the square of its argument if the 

synchronism condition is only approximately met. For this case, we 
show directly numerical results instead of using the approximate ana- 
j lysis in [19]. 

i 

In Figs. 5(a) - (d), we illustrate the numerical results for a three- 
coupled waveguide used as a power transfer device. The waveguide widths are 

i 

d^ * d 2 a d^ * 2 ym, the refractive indices are n^ * n 2 * n^ 53 3.44 , and the 
outside refractive index ng 33 3.436 [19]. The wavelength \ is 1.06 ym. The 
waveguides 1 edge-to-edge separations t are varied between 0.9 ym (near cutoff) 
to 4 # 4 ym. In Fig. 5(a), the propagation constants of the three supermodes using 



20 


the strongly coupled-mode theory are plotted and compared with those calculated 

exactly from solving the multilayered (slab) structure numerically. One finds 

very good agreement. A small discrepancy occurs for y^ w hen that mode is close 

to cutoff near t * 0.9 pm. The overlap integrals (solid curve) and 

(dashed curve) are plotted in Fig. 5(b) where aa large as 0.49, i.e., 

coupling is indeed very strong. (C^ 25 0.125 at t * 0.9 pm is also large). The 

coupling length L c - 2 tt/(y 1 “ Y 3 ) is plotted in Fig. 5(c) versus the waveguide 

separation t. One finds the extinction ratio p Qut ^ due t0 cross talks as shown 

in Fig. 5(d) decays as the waveguide separation t is increased. This has been 

discussed in [19] for a fixed separation t =* 1.9 pm using a different approach. 

Our result, at that separation, gives p Qut ^ ® 0*1082 = -9.66 dB which is 

actually higher than -12 dB given in [19] where the overlap integral between the 

two outer guides (= 0.0435) has been ignored. The results here should be 

more accurate since the exact propagation constants y , and at t * 1.9 pm 

are calculated accurately up to 7 digits after the decimal point and are also 

confirmed by the strongly coupled-mode theory. Taking the ratio of the extinc- 

2 

tion ratio P . to the square of the overlap integral, C 1#> , one finds 
out,l 12 

7 

P ,/C 10 « 0.78 at t 3 0.9 pm; 1.08 at t - 1.9 pm; 1.42 at t =* 2.9 pm, and 
out , 1 12 

1.75 at t = 3.9 pm. Thus one may only say that the extinction ratio is 

roughly proportional to the square of the overlap integral The propor- 

2 

tional constant estimated in [19] is * it /2 =* 4.9 for weak coupling and 

4.9/3 = 1.63 for strong coupling. The latter seems to agree better with our 

2 

results since the coupling is pretty strong here. Thus the factor tt~/ 2 is not 
appropriate for the example presented in [19]. The strongly coupled-mode theory 
should be applied when numerical accuracy is essential. 
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5. OPTICAL WAVELENGTH FILTERS 

Optical wavelength filters using waveguide couplers have been reported for 
Ti:LiNbO^ and InGaAsP - InP materials. The In^_ x Ga x ASyP^_y - InP material 
system is especially interesting because of its applications at 1.3 ym or 
1.55 pm wavelength and its potential for optoelectronic integrated circuits. 

The experiment reported in [8] has two coupled waveguides: one has a narrower 

guide width d - 0.42 ym, but a larger refractive index n (obtained following 
[21]) with y a ■ 0*127; the other has the guide width d^ * 0*91 ym and y^ - 0.078 
(n^ is also obtained following [21]). The Ga mole fraction x a (or x^) depends 
on the As mole fraction y a (or y fe ) for lattice matching [21]. The input power 
is assumed to be 1 in waveguide b. The results of a direct numerical approach 
have been shown in [8] and compared with the experimental data. We have applied 
the strongly coupled-mode theory using Eqs. (25) and (26) (exchanging a and b 
since the input is in guide b instead of in a) and compared our theoretical results 
with the experimental results in Fig. 6. The agreement is very similar to that 
in [8]. The parameters reported in [8] used for the theoretical calculations 
are within the measurement accuracy. No detailed explanations are given for the 
small discrepancy between the results for the theory and the experiment. We 
think the possible reasons may be (1) there is still some difference between the 
theoretical model in [21] and the experimental values for the refractive index, 
and (2) the losses in the waveguides are not taken into account. However, the 
comparison shown in Fig. 6 does show very good results. 
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6. CONCLUSIONS 

A strongly coupled-mode theory [9] - [11], [13], [14] has been presented and 
combined with the theory of Marcatili [15] for the two coupled waveguide case. 

The applications to two- and three-waveguide couplers, including power transfer 
devices and power dividers, have been investigated. This coupled-mode theory is 
applicable to very general cases for parallel dielectric waveguides with strong 
coupling and modes of general polarizations. It also accounts for the asymmetry 
of the waveguides and satisfies the energy conservation law and the reciprocity 
theorem self-consistently [13], [14]. The cross-talk problems in two and three 
coupled waveguides and their applications as optical wavelength filters have 
also been investigated and compared with the experimental data [8]. 
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APPENDIX A 

THE MATRIX ELEMENTS C , <T , K , T AND THE FIELD EXPRESSIONS 

pq pq pq pq 


(a) Matrix elements for the overlap integrals. 

C -~ // * H ( ^«zdxdy 

pq 2 Ji t t J 


C - 4- (C + C ) 
pq 2 pq qp' 


qp 


(b) Matrix elements for the coupling coefficients. 
"Conventional" coupling coefficients: 




ae <<l, - - e (,) (*,y) 


dxdy 


where s(x,y) is the permittivity function of the multiwaveguide system and 
e^(x,y) is the permittivity function of a single waveguide q. 

New Kpq used in Eq. (4) of this paper [13), [14]: 


k * t if 

pq 4 J [_t t 


e 22 


dxdy 


(c) The field expressions for the supennode. 
The transverse comoonents are 


e c - l a p( z > E t P>(x *y) 


(Al) 

(A2) 


(A3) 

(A4) 


(A5) 


(A6) 
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H - l a (z)H^ p) (x,y) . 

P P 

The longitudinal components are 

V^ P (2) r' E i P)(x,y ) 

P 

H z - l a p (z)H^ p) (x,y) . 


(A7) 


(A8) 

(A9) 
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APPENDIX B 

A FORMAL TREATMENT OF TWO COUPLED WAVEGUIDES FOR A LOSSLESS SYSTEM 
The coupled-mode equations are assumed to be of the form in general: 


~ a(z) - lY a a(z) + lK ab b(z) (Bl) 

^ b(z) - iY b b(z) + IK^aCz) (B2) 

where we have assumed for Che transverse fields 

B t - a(z)E^ a) + b(z)E^ b) (B3) 

H - a(z)H^ a) + b(z)H^ b) , (B4) 

t t t 


and the transverse components E^ a \ E^ b \ and H^ b ^ are all real. Thus, one 

finds that Y , Y. > K , and K, are all real. Power conservation leads to 
a D ab ba 

0 ‘h lM "3?T r ' II \ ” V 5d * d >- 


d « r * *•» *3* *1 

* dz Re aa + ab *21 + ba *12 * bb 1 

« 4— [aa* + (ab* + ba*)c“ + bb*] 
dz 

- .i*l<V - V + h»- *.b> 

* b ** I< V - V * K «b - K b.> 

where we have used Eqs. (Bl), (B2), the fact that C^» C 21 are real, and 


(B5) 


- C 12 + C 21 
C * o 


'12 


(B6) 
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Since both a and b are arbitrary, we conclude that the coefficients in front 
* * 

of ab and ba are zero and obtain: 


ab 


-S. 


< T . ' V c • 


(B7) 


which is the general lossless condition that the four parameters in the coupled- 
mode equations (Bl) and (B2) must satisfy. 
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FIGURE CAPTIONS 

Figure 1. (a) The asynchronism 5 is plotted versus the refractive index 

difference of two coupled waveguides, A n 3 n a - n^, which is 
proportional to the applied voltage V. 

(b) The overlap integrals C{2 (dashed line), C 21 (dotted line) and 
C * (C 12 + C2i)/2 (solid line) are shown. 

(c) The coupling coefficients K a ^, K^ a and /K a bK a b are plotted 
versus An. 

(d) The output powers in guide a, P a (solid curve), and in guide b, 
Pfc (dashed curve), are plotted versus An. The parameters are 
d a 3 d' D * 2 ym, waveguide edge-to-edge separation t * 1.9 um, 
wavelength X 3 1.06 ym. n a 3 2.2 + An/2, n^ 3 2.2 - An/2. 

The outside refractive index no 38 2.19. The coupler length 
£ 3 0.5811 mm. 

Figure 2. (a) The asynchronism 2^2 * Y1 ~ Y3 of three coupled waveguides are 

plotted versus the refractive index difference An 3 n^ - n2 
which is proportional to the applied voltage. The parameters 
are all similar to those in Fig. 1. di * d2 * d 3 3 2 ym. 
t * 1.9 ym, X 3 1.06 ym, n^ 3 n 3 3 2.2 + An/2, n2 3 2.2 - An/2. 
Loo 81 0.8105 mm. 

(b) The output powers PquT,1 (solid curve) and PoUT,3 (dashed curve) 
are shown for the power transfer devices with input power 

?IN * 1 in waveguide 1. 

(c) The output powers PquT,1 * P 0UT,3 are plotted versus An when the 
three waveguide couplers in (a) are used as a power divider. 

Figure 3. (a) The output powers PoUT,i (solid curve) and PoUT,3 (dashed curve) 

are plotted versus the coupling distance £ normalized to for 
a power transfer device. The same parameters from Fig. 2 are 
used except that An is fixed to be zero and £ is varying. 

(b) The output powers PoUT,l ** ^0UT,3 are plotted versus the nor- 
malized distance I/Lq when the three coupled waveguides in 
Fig. 2(c) are used as power dividers. (An 3 0, £ is varying 
here. ) 

Figure 4. (a) The output powers P()UT,1 (solid curve) and PquT, 3 (dashed curve) 

are plotted versus the normalized coupling distance L/Lq for the 
case An 3 -0.00023 where 2y2 ~ Y1 * Y3 * 0* 

(b) Similar conditions hold as 4(a) except that the device is used 
as a power divider. 
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Figure 5 


Figure 6 


(a) The propagation constants of the three supermodes using the 
strongly coupled-mode theory (dotted lines) are compared with 
the exact numerical calculations (solid lines). The waveguide 
edge-to-edge separation t is varied. The parameters are 

di * d2 * d3 * 2 yin, nj * n2 31 n^ = 3.44, nQ s 3.436. 

X * 1*06 ym. 

(b) The overlap integrals Ci 2 (*C 2 i) and Ci3(~C3i) are plotted versus 
the wavelength edge-to-edge separation t. 

(c) The coupling distance L q m 2tt/(y1 * Y3^ I 3 illustrated. 

(d) The output power PoUT 3 (hashed line) in the guide 3 and the 
extinction ratio Pqut]i (solid line) due to cross talk are 
shown. 

(a) The output powers of two coupled In x Gai- x AsyPi-y - InP wave- 
guides used as an optical wavelength filter. Tne input power is 
assumed to be 1 in waveguide b. The theoretical results for 
output powers at waveguide a, P a (solid curve), at waveguide b, 
Pb (dashed curve), are compared with the experimental data 
(circles) for P a . 
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ABSTRACT A strongly coupled-mode theory for reciprocal anisotropic 
multiwaveguide system is derived. The general reciprocal anisotropic 
medium is described by a symmetric permittivity tensor that can have non- 
zero off-diagonal elements. The derivation is based on the generalized 
reciprocity relation. The coupled mode equations are applicable to both 
lossy (gain) and lossless systems. For the special case of lossless 
systems, it is shown that the matrices in the coupled mode equation are 
Hermitian so that energy conservation is obeyed exactly. For the special 
case of a single anisotropic waveguide, our results also reduce to the pre- 
viously derived solutions by Marcuse. The strongly coupled-mode theory in 
an anisotropic multiwaveguide system is also illustrated with numerical 
examples. 
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1 . Introduction 

The application of the coupled mode theory to integrated optical devices 
such as waveguide couplers [1] — [3] , laser arrays [4] -[5] and optical filters 
[6]-[8] has been well known. Although it has long been recognized that 
the previous coupled mode theory is only applicable to weakly coupled 
systems [3]-[4], significant improvements for strongly coupled waveguides 
have only been developed recently in a series of papers £9]— [17] • The 
major improvement is the inclusion of the overlap integral C between 
modes of different waveguides that have been neglected in the conventional 
theory. 

Many of the integrated optical devices are made of materials that are 
anisotropic [18]-[22] . In anisotropic coupled waveguides, numerical solu- 
tions are usually difficult to calculate. In this paper we generalize the 
strongly coupled mode theory to an anisotropic medium in multiwaveguide 
systems. The strongly coupled mode equations are derived for a general 
reciprocal anisotropic medium. The general reciprocal anisotropic medium 
is described by a symmetric permittivity tensor that can have nonzero off- 
diagonal elements. The derivation is based on the generalized reciprocity 
relation [14]— [15] . In Section 2, the generalized reciprocity relation is 
extended to an anisotropic reciprocal medium. In Section 3, we briefly 
discuss the mode solutions and mode orthogonality relations for a general 
reciprocal anisotropic medium. In Section 4, the coupled mode equations for 
a general reciprocal anisotropic medium are derived for multiwaveguide 
systems. Expressions for the overlap integrals and the coupling coef- 
ficients are given. The derived coupled-mode equations are applicable to 
both lossy (gain) and lossless systems. In matrix notation, the coupled 
mode equations are of the form 


3 


cf 

dz 


i Q a 


where C is the overlap integral matrix. The elements of the vector a 
represent the amplitudes of the modes a^(z) in the waveguides, which can be 
modes of different polarizations in the same waveguides and/or different 
waveguides. Two special cases are examined next. In Section 5, the spe- 
cial case of a lossless system is studied. It is shown that in this case 

m a 

both C and Q are Hermitian matrices. Hence energy conservation is obeyed 
self-consistently. In Section 6, the special case of coupling of modes in 
a single waveguide is examined. In this case, the matrix C is diagonal. It 
is then shown that in this case the coupling coefficients are Identical to 
those previously derived by Marcuse [19]. In Section 7, the coupled mode 
theory is illustrated with a numerical example using Ti-diffused LiNbO^ 
waveguides. The application of voltages introduce: nonzero off-diagonal 
elements in the permittivity tensor. Thus there is coupling between the TE 
and TM modes in addition to coupling between modes of different waveguides. 
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2. Generalized Reciprocity Relation for Anisotropic Medium 


Let be the solution to Maxwell’s equations and the boundary 

conditions in a medium (x,y) for the whole space and E , and 

< 2 > 

H be the solutions to Maxwell’s equations and the boundary conditions in 
another medium (x,y). Both media are reciprocal so that the permit- 

tivity tensors satisfy the relations 


a<1>T r \ 
e (x,y) 

-<2>T, . 

e (x,y) 


; <l> (x.y) 

e <2> <x,y) 


where superscript T denotes transpose. 


a) 

( 2 ) 


It is straightforward to show that 


_ ,_<1> „<2> _<2> n <lX . _<1> ,=<2> »<1X _<2> 

V • (E x H -E xH )=i(i)E • (e -e ) • E 


(3) 


with the same procedure used for deriving the Lorentz reciprocity relation. 
The time convention exp(- iut) is used in this paper. When applied to a 
cylindrical geometry, the generalized reciprocity equation (3) becomes [14] 


3 _ 

3z 


n . , „ ,„<1> _<2> <2> _<1X 

/ J dxdy z • (E xfl -E xH ) 

, rr J J _< 1 > r< 2 >, , *< 1 >. 0 _< 2 > 

* im J / dxdy E • (e (x,y) - e (x,y)j • E 


( 4 ) 



3. Mode Solutions in Reciprocal Anisotropic Medium 
and Mode Orthogonality 
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In this section, we review the mode solutions and orthogonal relations 
in an anisotropic medium. In a reciprocal anisotropic medium, 


v 


XX 


xy 


xz 


'xy 


yy 


yz 


€ xz\ 


yz 


zz 


(5) 


/ 


The following decomposition is convenient for analyzing waveguide modes 


e =■ e + e,. + e+e 

t tz zt zz 

where subscript t denotes transverse components. 


( 6 ) 


e t * e xx ** + e xy + + e yy 99 


'tz 


e ftz + e £z 
xz yz 7 


e zt " e tz 


e • zz e 
zz zz 


(7) 

( 8 ) 

(9) 

( 10 ) 


Both e.. and e are symmetric while the transpose of e Is s and vice 

t ZZ CZ Z L 

versa. The longitudinal field components can be expressed in terms of the 


transverse field components as follows 
1 


H 


, V x E 

iu)U t t 


E 


a) e 


zz 


( 7 t * H t + iu «zt * M 


( 11 ) 


If e (x,y) is translational invariant in the z direction, then mode solu- 
tions exist. For modes propagating in the +z direction, the modal solutions 
of the electromagnetic fields are 
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E (q) (x,y) e i8< 1 z 
H (q) (x,y) e lB< l Z 

where q is the mode index. There also exists a mode with the same propa- 
gation constant and propagating in the -z direction. These are denoted 
by 

E ( " q) (x,y) e- iS «l z 
H (_q) (x,y) e' if3< l Z . 

For lossless systems, the two sets of mode functions (E^ q \ H^ q ^) and 

(E^ q \ q ^) have a simple relation (Section 5, equations (31)-(34)). 

However, for a general lossy (gain) reciprocal anisotropic medium, there is 

no simple relation relating the two sets of solutions [19], [23], (If 

e * 0 and e * 0, a simple relation such as that in [15] for z-inversion 
tz zt 

symmetry also exists,) 

Using the generalized reciprocity relation (4), mode orthogonality 
relations can be derived readily. Let 9 e» (E^\ ^) * 

(E (q) , H (q) ) exp (i0 q z), and (E <2> , H <2> ) - (E (p) , H (p) ) exp (i0 p z). It 
then follows that 

// dxdy 2 . (s£ q) x H< p) - E< p) x H< q) ) - 0 (12) 

for S * ~ 8 • Similarly, the following orthogonality relation 

p q 

If dxdy z . (E^ q) x H*" p) - b£" p) x H^ q) ) - 0 (13) 


also holds for 8 * 8 . 

q p 
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4. Coupled Mode Equations 


Consider a multiwaveguide system (Figure 1) with permittivity e (x,y), 
Also let e^ P \x,y) be the permittivity for a single waveguide p. We use 
the generalized reciprocity relation in the following manner. Let 

■<2> 


(x,y> - 

; (p) (*,y) 

(14) 

e <2> - 

(E (_p) + E ( " p) ) e" iS P z 
t z 

(15) 

h <2> - 

H < ' p> ) e" 18 P 2 
t z 

(16) 


be the permittivity and mode solution for a single waveguide p. For <1> we 
let 

e <l> (x,y) » e (x,y) (17) 

of the multiwaveguide system (Figure 1). The transverse field components 
of the multiwaveguide system are then approximated by linear superpositions 
of the individual waveguide modes of single waveguides. 

N 

q-l ^ 

N 


* 'i a (z) E^ q \x,y) 
t q t 


(18) 


“ l a (z) H^'Cx.y) 
C q=l q C 


.(q)< 


(19) 


where N is the number of modes that enter into the approximation. Note 
that each single waveguide q can contribute more than one mode (e.g., 

TEq and TMq modes). 

It then follows from Maxwell's equations that the longitudinal com- 
ponents are then given by the expansion (Appendix A) 

E <i> , l ,«> + <s<<> - : > . ,<«)) 

z L . e 1 zz z zt zt t 1 

q=»i zz 


( 20 ) 
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H a> - 

z 


l a fl (z) H 

q=*i 


(q) 


( 21 ) 


Substitute (14) - (21) into both sides of the generalized reciprocity rela- 
tion (4). The left-hand side of (4) gives 


N 

- I 

q-1 


- 4 - (a (z) e iS p Z ) 4 C 

dz q pq 


( 22 ) 


where 


C " J // dxdy z • [~E^(x,y) x P \x,y) 
Pq — « 

+ E^” p) (x,y) x H< q) (x,y)] 


(23) 


is the overlap integral which generally is not zero. The formulation with 
the neglect of the C terms with p * q has been referred to as the conven- 

pq 

tional coupled mode theory. The right-hand side of equation (4) gives 


N 


-i f 4 K a (z) e" iB P Z 

4* nn n 


q-1 


qp q 


(24) 


where 

K 

qp 


(e' p> - s t ) • «< q) 


+ E 


(-P). 


(e 


, , e(q) 

(p) = \ zz 


tz 


e tz> e 


zz 


E 


(q) 


(-p). 

( *t5 > - *«> 

• - :«> 

. E (q) 
t 

t 

e 

zz 

(-P). 

z 

- s«> 

-(q) 

* E t 


(-P). 

<4 P) - « ) 

ZZ 6 zz ; 

=(q) 

(q) 


zz 


zt zt 


+ c (-p) ^ £ zz^ e zz^ £ (q) E (q) , 
** * zz z * 


ZZ 


(25) 
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j Equating 


! 

! 


( 22 ) 

N 

l 

q-i 


and 

C 

pq 


(24) gives the coupled mode equation 
da (z) N 


S_ 

dz 


i T Q a (z) 
q-1 


I 

where 

i 

i 



K + 8 C 

qp p pq 


(26) 

(27) 


I 

The coupled mode equation (26) can be conveniently cast in matrix form, 

I 

C - i q a ( 2 ) , (28) 

> 

I 


where the matrices C and Q are of dimension N x N and the column vector a 
is of dimension N x 1, 
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5. Coupled Mode Equations for Lossless Multiwaveguide Systems 


The lossless (gain) medium is a special case of the general lossy 
system when the loss approaches zero. In this section, we shall derive the 
coupled mode equations for lossless systems. It will be shown that in this 
special case the C and Q matrices of the last section are both Hermitian so 
that energy conservation is obeyed exactly. 


For a lossless reciprocal anisotropic medium for both e and e 


(q) 


we 


have 


Let 


»T* 

e * e 

■(q)T* , |(q) 


(E (q) + E (q) ) e iS 9 z 


(H (q ) + H (q) ) e i0 9 z 

t Z 


(29) 

(30) 


be Che field solutions of a waveguide in a single waveguide 
Then by Caking complex conjugates of Che Maxwell equations, 

(E^* + E^* ) e" i8 9 z 
t z 

(-H (q) * - H (q) * ) e" 1 ^ 2 
t z 

is also a set of modal solutions and is propagating in the -z direction 
with a propagation constant g^. In this case, the modal fields propagating 
in the -z direction are simply related to those of the +z direction as 


q with real g . 

q 

it follows that 


follows 
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- „<«>* 


- -CO* 


Using (31) - (32) in (23), it follows that 


C pq m i If dxd y 1 


[g(q) x „<?)*+ E vp," x H vq ;j 


'n'l* r* - * 


Hence 


C - C 

pq qp 


and C is Hermitian 

. 

• 

j To show that 

Q is Hermitian, we make 

generalized reciprocity relation (4). Let 

=<1> 

-(q) 

e “ 

i 

e 

! E <1> . 

(E^ q ^ + E^) e i0< * z 
t z 

H <1> - 

I, 

(H (q) + H (q) ) e ifi q 2 

t Z 

and 

=<2> 

-(p) 

e * 

e 

e <2> - 

(E (p) *+ E (P) *) e" 10 P Z 
t z 

h <2> = 

(-H (p) *~ H (p) *) e‘ l8 P z 
t z 


( 42 ) 
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Then, from (4) It follows that 

jo 
4 


<*, - v ■ * t a d * d > E<p) * • 


- f!l**y (< p) *- (^ p) - J' q) ) • «< q) 


+ 

E (p) *. 

(-(p> 

- e (q) ) 

. E (q) 


t 

tz 

tz 

z 

+ 

e (p>\ 

<£ <p) 

- : <q> ) 

. R< q > 


z 

zt 

zt 

t 

+ 

E (p>* 

z 

(e (p) 

zz 

- 

E (q) 

z 


(43) 


Using (31) - (34) in (25), the quantity K - K can be calculated. 

qp pq 

After a moderate amount of algebraic manipulations, it can be shown that 
* 

K - K is equal to the right-hand side of (43). Hence, 

qp pq 


(3 - 8 ) C » K -K 

q p pq qp pq 


(44) 


It then follows from (27), (36) and (44) that 


Q - Q * - K +8C - K * - 8 C* 

pq qp qp p pq pq q qp 

- 0 (45) 

Hence, Q is Hermitian. It is straightforward to show that these relations 
(44) and (45) reduce identically to those of the isotropic waveguide case 
in [15). 

Energy conservation can be demonstrated as follows. The power P(z) in 
the raultiwaveguide system is 

P(z) - Y R e If dxdy z - E t x H* (46) 



Using (18) and (19) in (46) gives 


P(z) 


N 


l 

p>q-i 


* 


a a 


P q 


C 

pq 


(47) 


Using (28) and the fact that both C and Q are Hermitian, it follows that 
dP/dz ” 0 and the power P(z) is independent of z. Energy conservation is 
obeyed. 
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6 • Coupled Mode Equations for Single Waveguide System 

*(p) * ( q ) » 

In a single waveguide system, we have e ** e and e is the per- 
turbed permittivity that causes coupling among the modes in the single 
waveguide. The coupled mode equations for a single waveguide with a general 
reciprocal anisotropic medium were derived previously by Marcuse [19]. In 
this section, we shall show that for this special case, our results are 
identical to those in [19]. 

Let 



denote the unperturbed permittivity of 
E^^ are the mode solutions of e* with 
Hence they obey the mode orthogonality 
orthogonal relation, it is easy to see 
matrix. Hence 


(48) 

( p ) 

the single waveguide. Then E and 
propagation constants 0^ and 0^. 
relation of (13). Using the mode 
that C as given by (23) is a diagonal 


C * C 5 

pq P pq 


(49) 


where 5 is the Kronecker delta. The coupled mode equations become 

pq 

da 


C (— 2 - is a (z)) = i 7 K a (z) 

p v dz p p 1 ‘ qp 1 


(50) 


for p ■ 1, ..., N. Using (48) in (25), it follows that some of the terms 
are cancelled and the expression for K becomes 

qp 


15 


qp 


~f If dxd y { E t' p)# ( e t “ « t > • E t 


(q) 


+ k (-p). . <r - 


ZZ 


<e zt • St* • E 


<q> 


+ Ik J . E 

6 


zz 


'tz tz 


(q) 

z 


+ B <-P >.Ik ( -\ . 


(e - e ) • B (q) 
e,„ zt zt t 
zz 


♦ E (-p) Ik (e ’ . e ) B (q) } 

t e zz zz z 1 
zz 


(51) 


It can be shown (Appendix B) that (51) is identical to equation (46) 


of Marcuse [19] . 
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7 . Numerical Example 

In this section we illustrate the coupled mode theory developed in the 
previous sections with a numerical example using Ti-diffused LiNbO^ wave- 
guides. Consider two coupled slab waveguides, a and b, of thicknesses, 
d^ and d fe , and edge to edge separation, t. Voltages and are applied 
respectively to the two slabs. The uniaxial axis is assumed to be in the y 
direction for both guiding and substrate regions (Figure 2). In the 
following analysis we shall include four modes coupling (N»4)., The four 
modes are TE^, TM^ of guide a and TE^ and TMq of guide b, which are denoted 
as TEq , TM*, TEq and TMq with the corresponding amplitudes a^Cz) , a 2 (z), 
a^(z) and a ^(z) in the same order. Then we make the following substitu- 
tions. Let represent the permittivity profile of a single waveguide a and 
represent the permittivity profile of single waveguide b with zero 
voltages. 

for |x| < —■ 

for |x| > y ~ 

(52) 


r 




for guide a, where 


^ 2 
n 


oa 


0 

0 


ea 


0 

0 


oa 


0 

0 


0 


0 


2 

n 

o 


( 53 ) 
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For guide b 


r 


ob 


0 

0 


eb 


ob 


d h 

for |x - x b | < y ~ 


d a d h 

with x. ■ t + - — + — 
b 2 2 


for |x - x fe | > ~ 


(54) 


Then, let e 
■( 2 ) 


( 1 ) 


( 1 ) „( 1 ) 


e and (S' H 


) ■ (E, H) of the TE_ mode of e . Also 

0 a 


let e'"' - e and (E (2) , H (2) ) 

a 


(E, B) of the TM rt mode of e . Note that 
u a 


both ^ and 
(E (3 \ H (3 >) - 


=* ( 2 ) ■ =3 ( 3 ) 

e are equal to e . Similarly, e 


(E, H) of the TEq mode of and (E 


(4) 


. ;<*> 

H (4 >) 


and let 
(E, H) of 


the TMq mode of e^. 


The permittivity e is then taken to be the permittivity profile of the 

multiwaveguide system with voltages V and V. . Based on the charac- 

a d 

teristics of LiNbO^, Che permittivity e(x,y) can be approximated [22] by the 
following expressions (55) and (56). 
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Let 


2 4 y _y 

i - n r' -r 1 - 
oy oy 22 d 


Y 2 2 _1 

-r i. n n 

51 oy eY d 


, 9 v 

Y 2 2 _Y 

r’ n n 
51 oy ey d 


1 


eY 


o / v 

n 2 + n 4 r Y —I 

n • n o j 

oy oy 22 d^ 


with y ■ a, b. Then 


(55) 


^ ,» 
e_ 


s 


for x < 


for |x - x b | < — 


e c otherwise (56) 

In the numerical results to be illustrated, we have assumed the input 
wave to be the TM^ mode of guide a. Hence the coupled mode equations are 
solved and the initial conditions are set by letting a (z-O)*!, for q*2 and 

q 

0 for q-1,3,4. The mode functions are normalized such that C =1, 

qq 

q =* 1,2, 3, 4. In the calculations, we have tested the power conservation with 

equation (47). Power conservation is checked numerically to within 10 ^ 

percent. The reason for using high accuracy is because the cross- 

—4 

polarization powers can be very small (of the order 10 or less). Thus 
this high accuracy based on double precision is only for self-consistent 
check. The transfer efficiency of the two modulators of Figure 3 are con- 
sidered. Note that the two directional couplers have identical input 
guides but different output guides. The coupling section is of length l . 

By matching the field solution to the output guides, it then readily 
follows that [13], [17] 
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for the output power of directional coupler 1. Similarly, 
powers of directional coupler 2 are 

* ^ ( j ,L °3j C k3 a j ( ° ■k Ct) ) 

‘ °4j C \<*» 

- y // dxdy 2 • (E^ q) x h< p) *) 


9 TEq 


»TMo 


where 


pq 


(57) 

(58) 


the output 


(59) 

(60) 
(61) 


In the numerical results of Figures 4 to 6 we have used the following 

parameters: d a » db * 1.6 um, t * 1.6 yin, X - 0.633 ym, n 0 * 2.281, 

a b -I? 

n e » 2.195, n oa - n 0 b * 2.286, n ea ■ n e b * 2.2, r 22 * r 22 * x 10 m / v > 

a b -12 TMo TMq 

and r 51 * r $1 * 28.0 x 10 m/V. In Figure 4, P a and Pb are plotted 

as functions of the coupling length l for V a * Vb » 0. There are no off- 

m TEo TEq 

diagonal elements for e(x,y) in this case and P a * ?b * 0, so that there 
is no coupling to the TE modes. We note that when the output power in 
waveguide b is maximum at around l m 700 yra, the output power in waveguide 
a is not zero. This cross-talk is due to the overlap of fields between 
modes in waveguide a and waveguide b [17]. The computed overlap integral 
C 2 ^ in this case is equal to 0.1753076 which is outside the region of vali- 
dity of the conventional coupled-mode theory. The computed coupling coef- 
ficient K 24 - K 42 - 0.2284878 x lo” 2 . 
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In Figures 5 and 6 the output powers are plotted as functions of the 

modulating voltage V with V *» V and V * -V and l * 700 pm. For the para- 

a d 

meters chosen and a voltage range of -20 volts < V < 20 volts, there is 

appreciable transfer of power between the two guides. The results of 

Figure 6 indicate that the transfer of power to the TE modes is of the order 

of -40 dB to -30 dB with larger than because of the larger 

b a 

overlap of the fields between the TM^ mode of guide a and the TE^ mode of 
guide b. 
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8. Conclusions 


The strongly coupled-mode theory for a reciprocal anisotropic medium in 
multiwaveguide systems has been derived in this paper# Both lossy (gain) 
and lossless cases are considered. In general, the TE and the TM modes of 
separate waveguides are coupled due to the anisotropy of the medium. For 
the lossless case, the energy conservation is shown to be satisfied in this 
coupled mode formulation. For the case of a single anisotropic waveguide, 
we show that our results reduce identically to those of Marcuse [19]. The 
numerical results for two coupled anisotropic waveguides with both TE^ and 
TMq mode couplings have been illustrated. 
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Appendix A : Derivation of the Longitudinal Field Components 

To derive the expansion of the longitudinal field components as given 
by (20) and (21), we note that 


.<!> 




ZZ 


ft * < l> + iu St 


*l l> ) 


(Al) 


Substitute the expansions (18) and (19) in (Al) and note that for the qth 
mode of waveguide q 


*t * H c q) 


- iw e 


(q) . E (q) - i« e (q) E (q) 
Zt t ZZ z 


(A2) 


One then obtains (20). To derive (21), note that 
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7 x B (q) - ia>u H (q) 
t t z 
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Substituting (18) in (A3) and using (A4) gives (21). 



Appendix B : 


Comparison of in (51) with the Result of Marcuse [19] 


In comparing (51) with equation (46) of Marcuse [19], it is to be 
noted that the permittivity tensor elements are defined differently. 


-(m) “(m) 

C. T 6 


(Bl) 


where superscript (m) denotes Marcuse* Also 


«(m) - - 

e ■ e + c 
t t tz 


■(«) 


£ + £ 
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(B3) 


With the above substitution in equation (46) of [19], it can be shown, 

after a moderate amount of algebra, that the coupling coefficient K is 

qp 

identical to expression (51). 



24 


References 


[1] H. Kogelnik, "Theory of dielectric waveguides," in Integrated Optics , 

T. Tamir, Ed., Second Edition, Ch. 2. New York': Springer-Verlag, 1979. 

[2] A. Yariv, "Coupled-mode theory for guided wave optics,” IEEE J. 

Quantum Electron. , vol. QE-9, pp. 919-933, 1973. 

[3] H. F. Taylor and A. Yariv, "Guided wave optics," Proc. IEEE , vol. 62, 
pp. 1044-1060, 1974. 

[4] E. Kapon, J. Katz and A. Yariv, "Supermode analysis of phase locked 
arrays of semiconductor lasers," Opt. Lett. , vol. 9, pp. 125-127, 

1984. 

[5] A. Hardy and W. Streifer, "Analysis of phase array diode lasers," Opt. 
Lett . , vol. 10, pp. 335-337, 1985. 

[6] R. C. Alferness and R. V. Schmidt, "Tunable optical waveguide direc- 
tional coupler filter," Appl. Phys. Lett., vol. 33, pp. 161-163, July 
1978. 

[7] R. C. Alferness and J. J. Veselka, "Tunable Ti: LiNb 03 waveguide 
filter for long-wavelength (X ■ 1.3 - 1.6 ym) multiplexing/ 
demultiplexing," in Tech. Dig. Conf. Lasers and Electro-optics , 

(Anaheim, CA) , pp. 230-231, 1984. 

[8] B. Broberg, B. S. Lindgren, M. Oberg and H. Jiang, "A novel integrated 
optics wavelength filter in In Ga As P - In P," IEEE J. Lightwave 
Technol. , vol. LT-4, pp. 196-203, 1986. 

[9] A. Hardy and W. Streifer, "Coupled mode theory of parallel waveguides," 
IEEE J. Lightwave Technol. , vol. LT-3, pp. 1135-1146, 1985. 

[10] A. Hardy and W. Streifer, "Coupled modes of multiwaveguide systems and 
phase arrays," IEEE J. Lightwave Technol. , vol. LT-4, pp. 90-99, 

1986. 

[11] A. Hardy and W. Streifer, "Coupled mode solutions of multiwaveguide 
systems," IEEE J. Quantum Electron. , vol. QE-22, pp. 528-534, 1986. 

[12] H. A. Haus, W. P. Huang, S. Kawakami and N. A. Whitaker, "Coupled-mode 
theory of optical waveguides," IEEE J. Lightwave Technol., vol. LT-4, 

1986. 

[13] E. Marcatili, "Improved coupled-mode equations for dielectric guides," 
IEEE J. Quantum Electron. , vol. QE-22, pp. 988-993, 1986. 

[14] S. L. Chuang, "A coupled-mode formulation by reciprocity and a 
variational principle," IEEE J. Lightwave Technol., vol. LT-5, pp. 5-15, 

1987. 


[15] S. L. Chuang, "A coupled-mode theory for multiwaveguide systems 
satisfying the reciprocity theorem and power conservation,” IEEE J. 
Lightwave Technol. , vol. LT-5, pp. 174-183, 1987. 

[16] W. Streifer, M. Osinski and A. Hardy, "Reformulation of the coupled 
mode theory of multiwaveguide systems," IEEE J. Lightwave Technol. , 
vol. LT-4 , 1986. 

[17] S. L. Chuang, "Application of the strongly coupled-mode theory to 
integrated optical devices," submitted to IEEE J. Quantum Electron . 

[18] E. 6. Spencer, P. V. Lenzo and A. A. Ballman, "Dielectric materials 
for electro-optic, elasto-optic and ultrasonic device applications," 
Proc. IEEE , vol. 55, pp. 2074-2108, 1967. 

[19] D. Marcuse, "Coupled-mode theory for anisotropic optical waveguides, 
Bell Sys. Tech. J. , vol. 54, pp. 985-995, 1975. 

[20] R. A. Steinberg and T. G. Glallorenzi, "Modal fields of anisotropic 
channel waveguides," J. Opt. Soc. Am. , vol. 67, pp. 523-532, 1977. 

[21] 0. G. Ramer, "Integrated optic electro-optic modulator electrode 
analysis," IEEE J. Quantum Electron. , vol. QE-18, pp. 386-392, 1982. 

[22] A. Yariv and P. Yeh, Optical Waves in Crystals New York: Wiley Inter 
national, 1984. 

[23] L. B. Felsen and N. Marcuritz, Radiation and Scattering of Waves 

Englewood Cliffs, New Jersey: Prentice-Hall, 1973. 


26 


Figure Captions 


Figure 1 (a) 

(b) 

(c) 


A single waveguide p described by £^ P \x,y) for the entire 
space 

A single waveguide q described by e^ q \x,y) for the entire 
space 

A multiwaveguide system described by e(x,y) in the entire 
space. 


Figure 2 


Geometric configuration of two coupled anisotropic wave- 
guides for the numerical example in Section 7. 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


Directional couplers with identical input guides a but dif- 
ferent output guides a (directional coupler 1) or b 
(directional coupler 2). 

The output powers in the TMq mode of guide a, P^^ » and in 
the TM rt mode of guide b, pP^, as functions of the coupling 

U D 

length l without applied voltage, V ■ 0. 

The output powers in the TM. mode of guide a, and in 

U 3 

the TM a mode of guide b, pP^, as functions of the applied 
U D 

voltage V » - V * V. 
a d 

The output powers in the TE mode of guide a, P^ E 0 , and in 

U 3 l 

TE mode of guide b, pP^, as functions of the applied 
U D 

voltage V = - V. = V. 
a b 
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